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ABSTRACT 

We investigate the launching of outflows from the disk-magnetosphere boundary of 

slowly and rapidly rotating magnetized stars using axisymmetric and exploratory 
4^ 3D magnetohydrodynamic (MHD) simulations. We find long-lasting outflows in both 

cases. (1) In the case of slowly rotating stars, a new type of outflow, a conical wind, is 
Q found and studied in simulations. The conical winds appear in cases where the mag- 

;-H netic flux of the star is bunched up by the disk into an X-t3^e configuration. The winds 

c/^ have the shape of a thin conical shell with a half-opening angle ^ 30° — 40°. About 

10 — 30% of the disk matter flows from the inner disk into the conical winds. The 

conical winds may be responsible for episodic as well as long-lasting outflows in dif- 
^ ferent types of stars. (2) In the case of rapidly rotating stars (the "propeller regime"), 

^ a two-component outflow is observed. One component is similar to the conical winds. 

A significant fraction of the disk matter may be ejected into the winds. A second com- 
Q^ ponent is a high-velocity, low-density magnetically dominated axial jet where matter 

CO flows along the opened polar field lines of the star. The jet has a mass flux about 10% 

CO that of the conical wind, but its energy flux (dominantly magnetic) can be larger than 

the energy flux of the conical wind. The jet's angular momentum flux (also dominantly 
Q magnetic) causes the star to spin-down rapidly. Propeller-driven outflows may be re- 

Q\ sponsible for the jets in protostars and for their rapid spin-down. The jet is collimated 

O by the magnetic force while the conical winds are only weakly collimated in the sim- 

ulation region. Exploratory 3D simulations show that conical winds are axisymmetric 
. i-H about the rotational axis (of the star and the disk), even when the dipole field of the 

^ star is significantly misaligned. 
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1 INTRODUCTION 

Outflows or jets are observed from many disk accreting ob- 
jects ranging from young stars to systems with white dwarfs, 
neutron stars, and black holes (e.g., Livio 1997). 

A large body of observations exists for outflows from 
young stars at different stages of their evolution, ranging 
from protostars, where powerful collimated outflows are ob- 
served, to classical T Tauri stars (CTTSs), where the outflows 
are weaker and often less collimated (see review by Ray et 
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al. 2007). Correlation between the disk and jet power had 
been found in many CTTSs (e.g., Cabrit et al. 1990; Har- 
tigan, Edwards & Gandhour 1995). A significant number of 
CTTSs show signs of outflows in spectral lines, in particular 
in He I where two distinct components of outflows had been 
found (Edwards et al. 2003, 2006; Kwan, Edwards, & Fis- 
cher 2007). Outflows are also observed from accreting com- 
pact stars such as accreting white dwarfs in symbiotic binaries 
(e.g., Sokoloski & Kenyon 2003), or from the vicinity of neu- 
tron stars, such as from Circinus X-1 (Heinz et al. 2007). 

Different theoretical models have been proposed to ex- 
plain the outflows from protostars and CTTSs (see review by 
Ferreira, Dougados, & Cabrit 2006). The models include those 
where the outflow originates from a radially distributed disk 
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Figure 1. Two-component outflows observed in slowly (left) and rapidly (right) rotating magnetized stars for the reference runs described in this 
paper. The background shows the poloidal matter flux Fm = pvp, the arrows are the poloidal velocity vectors, and the lines are sample magnetic 
field lines. The labels point to the main outflow components. 



v\Aind (Konigl & Pudritz 2000; Casse & Keppens 2004; Ferreira 
et al. 2006) or from the innermost region of the accretion 
disk (Lovelace, Berk & Contopoulos 1991). Further, there is 
the X-wind model (Shu et al. 1994; 2007; Najita & Shu 1994; 
Cai et al. 2008) where most of the outflow originates from 
the disk-magnetosphere boundary. The maximum velocities 
in the outflows are usually of the order of the Keplerian ve- 
locity of the inner region of the disk (or higher) . This favors 
the models where the outflows originate from the inner disk 
region, or from the disk-magnetosphere boundary (if the star 
has a dynamically important magnetic field) . 

Outflows from the disk-magnetosphere boundary were 
investigated in early simulations by Hayashi, Shibata & Mat- 
sumoto (1996) and Miller & Stone (1997). A one-time 
episode of outflows from the inner disk and inflation of the 
innermost field lines connecting the star and the disk were ob- 
served for a few dynamical time-scales. Somewhat longer sim- 
ulation runs were performed by Goodson et al. (1997, 1999), 
Hirose et al. (1997), Matt et al. (2002) and Kiiker, Henning 
& Riidiger (2003) where several episodes of field inflation 
and outflows were observed. These simulations hinted at a 
possible long-term nature for the outflows. However, the sim- 
ulations were not sufficiently long to establish the behavior 
of the outflows. Much longer simulation runs were obtained 
by treating the disk as a boundary condition (e.g. Fendt & Ei- 
sner 1999, 2000; Matsakos et al. 2008; Fendt 2009; see also 
von Rekowski & Brandenburg 2004; Yelenina, Ustyugova & 
Koldoba 2006). These simulations help understand, for exam- 
ple, the roles of the disk wind and stellar wind components 
in the outflow and collimation. However, for understanding 
the launching mechanisms it is important to have a realistic, 
low-temperature disk and to solve the full MHD equations in 
all of the disk and coronal space. 

The goal of this work is to obtain long-lasting (robust) 
outflows from a realistic low-temperature disk (not a bound- 
ary condition) into a high-temperature, low-density corona. 
We obtained such outflows in two main cases: (1) when the 
star rotates slowly but the field lines are bunched up into an 
X-type configuration, and (2) when the star rotates rapidly, 
in the propeller regime (e.g., Illarionov & Sunyaev 1975; Al- 
par & Shaham 1985; Lovelace, Romanova & Bisnovatyi-Kogan 
1999) and the condition for bunching is also satisfied. In 
both cases, two-component outflows have been observed (see 
Fig. [TJ. One component originates at the inner edge of the 
disk and has a narrow-shell conical shape, and therefore we 



call it a "conical wind". The other component is a magneti- 
cally (or centrifugally) driven high-velocity low-density wind 
which flows along stellar field lines. We call it a "jet". The jet 
may be very powerful in the propeller regime. Below we dis- 
cuss both regimes in detail (see §3 -§4) after description of 
the numerical approach (see §2). In §5 we discuss different 
properties of outflows. In §6 we present exploratory 3D simu- 
lations of conical winds, and in §7 we compare conical winds 
and propeller outflows with the X-wind model. In §8 we ap- 
ply the model to different types of stars, and in §9 we present 
our conclusions. Appendixes A and B clarify different aspects 
of the numerical model. Appendix C summarizes results of dif- 
ferent runs for a variety of parameters. 



2 NUMERICAL MODEL 

We simulate the outflows resulting from disk-magnetosphere 
interaction using the equations of axisymmetric MHD de- 
scribed below. Axisymmetric simulations of the outflows are 
similar to those performed earlier for the propeller regime 
(e.g., U06), but differ in initial and boundary conditions. Be- 
low we give an outline of the numerical model. 



2 . 1 Basic Equations 

Outside of the disk the flow is described by the equations of 
ideal MHD. Inside the disk the flow is described by the equa- 
tions of viscous, resistive MHD. In an inertial reference frame 
the equations are: 



dp 



dt 

dt 
OB 



+ V • (pv) = , 



dt 
dt 



+ v-r = pg. 



V X (v X B) + V X (rytV x B) = , 



+ V • {pS^r) = Q , 



(1) 



(2) 



(3) 



(4) 



Here, p is the density and S is the specific entropy; v is the 
flow velocity; B is the magnetic field; T is the momentum 
flux-density tensor; Q is the rate of change of entropy per 
unit volume; and g = —{GM/r'^)r is the gravitational accel- 
eration due to the star, which has mass M. The total mass of 
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Figure 2. Conical winds at different times T. The background shows the poloidal matter flux pvp (with the scale below the plots), the arrows 
are the poloidal velocity vectors, and the lines are sample magnetic field lines (same set in all frames). Time T is measured in Keplerian rotation 
periods at r = 1. For example, for CTTS To = 1.04 days (see Table 1) and time T = 700 corresponds to 2 years. The sample vector Vp = 1 
corresponds to vq = 195 km/s. 
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Figure 3. The poloidal matter flux pvp (with the scale on the right-hand side), sample magnetic field lines, and velocity vectors in the conical 
wind at time T = 500. Sample numbers are given for the dimensionless poloidal Vp and total vt velocities, and for the density p. To obtain 
the dimensional values one needs to multiply these numbers by the reference values given in the Table 1. For example, for application to CTTS: 
Vp = 1 corresponds to vq = 195 km/s, p = 1 corresponds to po = 4.1 x 10~^^ g cm~^, and the distance r = 1 corresponds to Rq = 0.02 AU. 



the disk is negligible compared to M. The plasma is consid- 
ered to be an ideal gas w^ith adiabatic index 7 = 5/3, and 
S — \n(p/p^). We use spherical coordinates (r.O.cj)) with 
measured from the symmetry axis. The condition for axisym- 
metry is d/dcj) = 0. The equations in spherical coordinates are 
given in U06. 

The stress tensor T and the treatment of viscosity and 
diffusivity are described in Appendix A. Briefly, both the vis- 
cosity and the magnetic diffusivity of the disk plasma are con- 
sidered to be due to turbulent fluctuations of the velocity and 



the magnetic field. We adopt the standard hypothesis w^here 
the molecular transport coefficients are replaced by turbulent 
coefficients. To estimate the values of these coefficients, we 
use the a-model of Shakura and Sunyaev (1973) where the 
coefficient of the turbulent kinematic viscosity z/t = avcl/Qx, 
where Cs is the isothermal sound speed and Qxir) is the Kep- 
lerian angular velocity. Similarly, the coefficient of the turbu- 
lent magnetic diffusivity rjt = adcl/QK- Here, ay and ad are 
dimensionless coefficients which are treated as parameters of 
the model. 
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Figure 5. The background shows different parameters of a conical wind at time T = 500, in the full simulation region, panels a-d, and near 
the star, panels e-h. The background shows (from left to right): the poloidal matter flux, Fm = pvp, the angular velocity Q, the poloidal current 
jp = rBcj), and the entropy S. Lines are sample magnetic field lines. Vectors are velocity vectors. The thick red line and the marks A, B, C show 
one of the field lines and positions along this line which are used for analysis of the forces and velocities in Fig. [6| The white vector in panel g 
shows schematically the direction of the magnetic force, M. 




Figure 4. Left panel: equatorial distribution of the poloidal magnetic 
field in the inner part of the simulation region at different times T. 
Right panel: same as left panel but for density. 



2.2 Reference Units 

The MHD equations are solved in dimensionless form so that 
the results can be readily applied to different accreting stars 
(see §7). We take the reference mass Mo to be the mass M 
of the star The reference radius is taken to be twice the ra- 
dius of the star, Ro = 2R^ . The surface magnetic field is 
different for different types of stars. The reference velocity is 
vo = [GM / RqY^'^ , The reference time-scale to = Ro/vo, and 
the reference angular velocity Qo = 1 /to . We measure time in 
units of Po — 27rto (which is the Keplerian rotation period at 
r — Ro). In the plots we use the dimensionless time T = t/Po. 
The reference magnetic field is Bo = B^{R^/Ro)^ /jl, where 
fl is the dimensionless magnetic moment. The reference den- 
sity is taken to be po = Bq/vq. The reference pressure is 
Po = Bq. The reference temperature is To = po /Tlpo — vl /IZ, 
where IZ is the gas constant. The reference accretion rate is 
Mo — povoRo' The reference energy flux is Eo — MoVo- The 
reference angular momentum flux is Lo — MovoRo- 

The reference units are defined in such a way that the 



dimensionless MHD equations have the same form as the di- 
mensional ones, equations (l)-(4) (for such dimensionaUza- 
tion we put GM — 1 and IZ — I). Table 1 shows exam- 
ples of reference variables for different stars. We solve the 
MHD equations (l)-(4) using normalized variables: p — p/ po, 
V = v/vo, B — B/Boy etc. Most of the plots show the normal- 
ized variables (with the tildes implicit) . To obtain dimensional 
values one needs to multiply values from the plots by the cor- 
responding reference values from Table 1 . 

2.3 Initial and Boundary Conditions 

We assume that the poloidal magnetic field of the star is an 
aligned dipole field B = [3(/l/ • r)r — \ir^\lr^ , where \i is 
the star's magnetic moment. The initial density and tempera- 
ture distributions are different in cases of conical winds and 
propeller-driven winds. 

23 A Conical Winds 

Initial Conditions. At time T = the simulation region is filled 
with a low-density, high-temperature, isothermal plasma re- 
ferred to as the corona. Initially it is non-rotating. Thus the 
density and pressure distributions in the corona are: 

p^pc exp[G'M/(7^Ter)], p ^ pc exp[G'M/(7^Ter)] , 

where Tc is the corona temperature, pc is the coronal density 
at the external boundary, pc = pcTlTc. 

We divide the external boundary Rout into a disk re- 
gion, Od < < and a corona region, < < Od, 
with 6d ~ 65°. Initially, there is no disk in the simulation re- 
gion. When the simulations start, we permit high-density low- 
temperature matter to enter the simulation region through 
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Table 1. Reference values for different types of stars. We choose the mass M, radius i?*, equatorial magnetic field B^ and the period P* of the 
star and derive the other reference values (see §2.2). To apply the simulation results to a particular star one needs to multiply the dimensionless 
values from the plots by the reference values from this table. 



the disk boundary region, > Od, with a fixed density p = pd- 
Matter continues to flow inward due to viscosity (see Ap- 
pendix A) . We increase the spin of the star gradually from 
a small value corresponding to rcor = Rout (where rcor = 
(GM/n^y^^) up to a final value Q*. Information about the 
stellar rotation propagates rapidly (at the Alfven speed) into 
the low-density corona. 

We did simulations for a variety of parameters. However, 
we take one case with typical parameters to be our refer- 
ence case and show the results for this case. In the reference 
case, the dipole moment of the star /i = 10; the density in 
the corona pc = 0.001, the density in the disk pd = 10; 
the corona is hot with temperature Tc = 1; the disk is cold 
with temperature Td = {pc/pd)Tc = 10~^. The angular ve- 
locity of the star corresponds to a corotation radius rcor — 3, 
= (GM/vcor) = 0-19- The coefficients of viscosity and dif- 
fusivity are ay = 0.3 and ad = 0.1. The dependences of our 
results on different parameters are discussed in Appendix C. 

The boundary conditions at the inner boundary r = Rin 
are the following: The frozen-in condition is applied to the 
poloidal component Bp of the field, such that Br is fixed 
while Be and obey "free" boundary conditions, dBe /dr = 
and dBcj, /dr = 0. The density, pressure, and entropy also 
have free boundary conditions, d{...)/dr = 0. The velocity 
components are calculated using free boundary conditions. 
Then, the velocity vector is adjusted to be parallel to the mag- 
netic field vector in the coordinate system rotating with a star. 
Matter always flows inward at the star's surface. Outflow to a 
stellar wind is not considered in this work. 

The boundary conditions at the external boundary r = 
Rout in the coronal region < 6> < 6>d are free for all hydro- 
dynamic variables. However, we prevent matter from flowing 
into the simulation region from this part of the boundary. We 
solve the transport equation for the flux function ^ so that the 
magnetic flux flows out of the region together with matter. If 
the matter has a tendency to flow back in, then we fix ^. In 
the disk region, Od < < 7r/2, we fix the density at /> = pd, 
and establish a slightly sub-Keplerian velocity, = /^O(rd), 
where k, = 1 — 0.003 so that matter flows into the simula- 



tion region through the boundary. The inflowing matter has a 
fixed magnetic flux which is very small because Rout ^ Rin- 
The boundary conditions on the equatorial plane and on 
the rotation axis are symmetric and antisymmetric. 

2.3.2 Propeller Regime 

The initial and boundary conditions for the propeller regime 
are the same as those used in R05 and U06. Here, we summa- 
rize these conditions. 

Initial Conditions. We place both the disk and the corona 
into the simulation region. We assume that the initial flow is 
barotropic with p — p{p), and that there is no pressure jump 
at the boundary between the disk and corona. Then the initial 
density distribution (in dimensionless units) is the following: 

{p/lZTd , p > Pb and rsinO > rt , 
p/lZTc , p < Pb or r sin < rb , 

where pb is the pressure on the surface which separates the 
cold matter of the disk from the hot matter of the corona. 
On this surface the density jumps from pb/Td to pb/Tc. Here 
Tb is the inner disk radius. Because the density distribution is 
barotropic, the angular velocity is constant on coaxial cylin- 
drical surfaces about the axis. Consequently, the pressure 
distribution may be determined from the Bernoulli equation, 

F{p) -\-^-\-^c = E = const . 

Here, $ = —GM/\r\ is gravitational potential, $c = 
/r'Tine* ^^(^)^^^ centrifugal potential, which depends only 
on the cylindrical radius r sin 0, and 

{7lTd\n{p/pb) , p> Pb and rsinO > Vb , 
IZTc \a{p/pb) , p <Pb or rsmO <rb . 

The angular velocity of the disk is slightly sub-Keplerian, 
VL{e = 7r/2) = kQk = I - 0.003), due to which the 
density and pressure decrease towards the periphery. Inside 
the cylinder r < Vb the matter rotates rigidly with angular 
velocity Q.(rb) = ^{GM/rlY^'^. For a gradual start-up we 
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Figure 6. Panel a shows the projection of forces onto part AB of the bold field line shown in Fig. [sj. The labels are: G— gravitational, 
C— centrifugal, M— magnetic, P— pressure gradient, and i?— sum of all forces. Panel b shows the projection of forces onto part BC of the field 
line. Panel c shows velocities along part AB of the field line: ^p— poloidal velocity, i;^— azimuthal velocity, ^t— total velocity, Alfven speed, 
'Usrn— slow-magnetosonic speed, v fm— ^SiSt- magnetosonic speed, and i;esc— escape velocity. Panel d shows velocities along part BC of the field 
hne. 




Figure 7. Panel a: Two components of outflow in the reference case [rcor = 3). Panel b: Same but for a more rapidly rotating star [rcor = 2). 
Both cases are shown at time T = 250. 



change the angular velocity of the star from its initial value 
Q{rb) = 5"^/^ ^ 0.09 in = 5) to a final value of 1^* = 1 over 
the course of three Keplerian rotation periods at r = 1. 

For the propeller regime we use a slightly different set 
of parameters compared w^ith the conical w^ind case (in order 
to be consistent w^ith our earUer simulations in R05, U06). 
Below^ w^e describe the similarities and differences: the dipole 
moment of the star is the same, /i = 10. The angular velocity 
of the star in the propeller regime is larger, = 1. The initial 
density in the disk (at the inner edge) is pd = 1, which is 
smaller than the external density ipd = 10) in conical winds. 
The initial temperatures are a factor of two smaller than in 
conical winds, Tc = {p/np)c = 0.5, Td = 0.0005. 

Boundary conditions for the propeller regime are similar to 
those for conical winds with the following differences. At the 
external boundary (disk region) we take free conditions for 
all variables. There is no condition of fixed density at the disk 
part of the boundary. 

The system of MHD equations ( [1|4| ) was integrated nu- 
merically using the Godunov-type numerical scheme (see Ap- 
pendix B). The simulations were done in the region Rin < 



r < Rout, < < n/2. The grid is uniform in the 0- 
direction. The size steps in the radial direction were chosen 
so that the poloidal-plane cells were curvilinear rectangles 
with approximately equal sides. A typical region for investi- 
gation of conical winds was 1 < r < 16, with grid resolution 
Nr X Ne = 51 X 31 cells. A typical region for investigation 
of the propeller regime was 1 < r < 48 with grid resolution 
Nr X Ne = 85 X 31 cells. Test simulations at angular grids 
Ne = 51 and A^^ = 71 were also performed. Each simulation 
run at the lowest resolution takes about two months of com- 
puting time on a single processor. We performed 40 different 
simulation runs for different parameters on our local cluster 
of 20 computers for the investigation of conical winds. 3D 
simulation runs were performed with "Cubed sphere" parallel 
code on NASA high-performance facilities. 
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Figure 8. Same as in Fig. [2] but for the propeller regime. The simulation time T = 2200 corresponds to 6 years. The sample vector Vp 
corresponds to = 2 x i;o = 390 km/s. 
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Figure 9. Same as in Fig.[3]but for the propeller regime at time T = 1400. 



3 MATTER FLOW IN CONICAL WINDS AND IN THE 
PROPELLER REGIME 

3.1 Matter flow, velocities, and forces in conical winds 

A large number of simulations were done in order to under- 
stand the origin and nature of conical winds. All of the key 
parameters were varied in order to ensure that there is no spe- 
cial dependence on any parameter (see Appendix C) . We ob- 
served that the formation of conical winds is a common phe- 
nomenon for a wide range of parameters. They are most per- 
sistent and strong in cases where the viscosity and diffusivity 
coefficients are not very small, > 0.03, > 0.03. Another 
important condition is that av > ad', that is, the magnetic 
Prandtl number of the turbulence, Prm = ay /ad > 1. This 
condition favors the bunching of the stellar magnetic field by 
the accretion flow 



For a discussion of the physics of conical winds we focus 
on one set of parameters which serves as a reference case. 
These parameters are: av = 0.3 and ad = 0.1; = 0.19 
(ircor = 3), /i = 10, pd = 10, and pc = 10-^ Tc = 1, 
Td — 10~^. The simulations were done in dimensionless form 
and can be appHed to different stars (see Table 1). How- 
ever, for illustration we often show dimensional examples for 
CTTSs with parameters taken from Table 1. For example, in 
application to CTTSs, r2* = 0.19 corresponds to P* = 5.4 
days and the unit of time used in the figures is Po = 1-04 days 
(see Table 1 for other reference values) . 

Fig. [2] shows snapshots of simulations at different times 
T. One can see that the cold dense disk matter enters the 
simulation region from the external boundary and moves in- 
ward towards the star on the viscous time-scale. The accre- 
tion flow bunches up the field lines of the dipole field to a 
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relatively small region near the star. All field lines shown at 
T = are bunched up close to the star by T = 100. The 
inclined configuration of the resulting poloidal field and in- 
flation of external field lines create conditions favorable for 
matter outflow from the inner disk. The outflow starts at 
T ^ 120 and gets stronger later. Matter flows from the in- 
ner disk into hollow, conical shaped winds with half-opening 
angle ~ 30° — 40°. The conical winds are non-stationary 
showing variations associated with events of inflation and 
reconnection of the magnetic field lines (see animations at 
http://www. astro. Cornell edu/^romanova/conical htm) . The 
simulation runs continue for a long time, about T = 740, 
which is about 2 years for CTTSs. The outflows remain strong 
until the end of the simulation runs. It is reasonable to con- 
clude that these accretion-driven outflows into the conical 
winds will persist as long as matter is supplied from the disk. 

Fig. [s] shows the configuration at T = 500. One can see 
that the disk matter comes close to the star and accretes onto 
the star through a small dense funnel. Some field lines are 
strongly inflated, and the conical wind flows from the disk 
along these lines. There is also a set of partially inflated field 
lines (a dead zone) where accretion does not occur (e.g., Os- 
triker & Shu 2005; Spruit & Taam 1990). Matter in conical 
winds rotates with the Keplerian velocity at the base of out- 
flow, Vcf) ~ vk- It continues to rotate rapidly in the conical 
wind at larger distances from the star. The poloidal velocity 
Vp increases gradually from very small values at the begin- 
ning of the outflow, up to Vp ^ O.bvK- The main contribution 
to the total velocity vt comes from the azimuthal component. 
There is another, high-velocity component of the low-density 
matter which flows along the stellar field lines. In application 
to CTTSs the velocity is > 200 km/s. 

Fig.[4|shows the variation of the density and the poloidal 
magnetic field along the equator in the inner part of the sim- 
ulation region at different times T. One can see that the disk 
matter has approximately constant density at different radii, 
but there is a density peak closer to the star. The peak in- 
creases with time, but it does not appreciably influence the 
fluxes calculated at the surface of the star (see Fig. [14]) . Fig. 
|4] also shows that the poloidal magnetic field of the star is 
compressed, and the compression increases with time. Com- 
pression of the star's magnetic field by the accretion flow is 
also assumed in the X-wind model (e.g., Najita & Shu 1994). 

To understand the physics of conical winds in greater de- 
tail we show in Fig.jsjthe distribution of different parameters 
at time T — 500. Panels h and / show that the innermost 
region of the closed magnetosphere rotates with the angu- 
lar velocity of the star (1.2 < r < 1.8). At larger distances 



(r > 2.5), the corona above the disk rotates with the an- 
gular velocity of the disk. Strongly inclined field lines which 
start in the disk go through regions of lower and lower an- 
gular velocity and are strongly wound up owing to the differ- 
ence in the angular rotation rates along the lines. This leads 
to a strong poloidal current flow Jp oc rB^ above the disk 
(see panels c and g) which gives rise to the magnetic force 
Fp ~ —V[{rB(pY]. This is the main force driving matter into 
the conical wind. Driving of winds by the magnetic force from 
the inner disk was proposed earUer by Lovelace et al. (1991). 
The direction of the magnetic force is shown schematically in 
Fig. [S^. It acts upwards and towards the axis, which is differ- 
ent from the centrifugal force. This determines three of the 
important properties of conical winds: (1) their small open- 
ing angle; (2) the fact that the wall of the cone is narrow, 
and (3). the gradual collimation of conical winds. If the cen- 
trifugal force were dominate (e.g. Blandford & Payne 1982) 
then the cone would have a wider opening angle and outflow 
would flow over a wide range of directions (as in the X-wind 
model of Shu et al. 1994). Panels d and h show the distribu- 
tion of entropy S which shows that matter flowing from the 
disk into the wind is cold and it is not thermally driven. To 
analyze the forces driving matter into the conical winds we 
select one of the field lines, s, (see red bold line in panels e- 
K) and we project forces onto this field line. We split the line 
into two parts (see panel a) . Part AB starts from the disk and 
ends at the place where the line curves towards the star; part 
BC continues from there to the surface of the star . 

Fig. |6]i shows the projection of all forces onto part AB 
of the field line. One can see that the main force accelerat- 
ing matter into the conical wind is the magnetic force M. The 
centrifugal (C) and gravitational {G) forces approximately 
compensate each other and the sum C + G is negative. The 
pressure gradient force P is small. The ^-component of the 
magnetic force leads to frequent forced reconnection events 
of the inflated field lines and to ejection of plasmoids into the 
conical wind. Panel h shows the projection of the forces onto 
segment BC of the field line. One can see that it is chiefly the 
centrifugal force which accelerates the low-density matter to 
high velocities in this region. Panel c shows that in the conical 
wind the poloidal velocity Vp (along part AB) gradually in- 
creases and crosses the slow magnetosonic {v — Vsm), Alfven 
{v — va) and fast magnetosonic {v — Vfm) surfaces. Mat- 
ter rotates rapidly, therefore the azimuthal component Vcf^ is 
much larger than poloidal one, and the total velocity vt is de- 
termined by the azimuthal rotation of the flow. Panel d shows 
that there is an interval of high velocity along the stellar part 
(BC) of the field line. Thus, we observe a two-component 
flow: (1) a high-density low- velocity conical wind which is the 
main component of the outflows, and (2) a low-density fast 
outflow along the stellar field lines which occupies a much 
smaller region. 

We find that the region of the fast coronal flow increases 
in size with the star's rotation rate. As an example we de- 
creased the corotation radius from Vcor — 3 (Q* = 0.19) to 
Tcor — 2 (^1* = 0.35) and observed that the region of fast 
coronal flow increased significantly. Fig. [t] shows the differ- 
ence. The region is even larger for smaller corotation radii 
when the star is closer to the propeller regime. In the pro- 
peller regime (see §4) the fast jet component occupies the en- 
tire region within the conical wind and is very powerful. The 
star spins up for both Rcor — 2 and 3. Cases Rcor — 1, 1.5 
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Figure 11. Forces (top panels) and velocities (bottom panels) along the field lines in the propeller regime atT = 1400. Left panels: forces and 
velocities along the closed magnetic field line which starts in the disk (at r = 4.3) and ends on the star. We take only the part of the line up to 
the neutral point where the line curves towards the star(Br = 0). Forces are projected onto the field line. Labels are the same as in Fig. [6] Right 
panels: forces and velocities along the open field Hne which starts on the surface of the star. Note that the scale for va, Vsm, 'f^fm is different from 
that for other velocities. 



correspond to the propeller regime where the star spins down 
due to the interaction with the disk and corona. We did not 
perform a refined search for the rotational equilibrium state, 
in which the star has alternate spin-up and spin-down peri- 
ods, but zero torque on an average (e.g. R02, Long et al. 
2005). In this state we expect the jet component to occupy 
a large part of the region above conical winds. Even a weak 
stellar wind (not considered in this paper) may enhance the 
jet component. 

3.2 Matter flow, velocities and forces in the propeller 
regime 

Here we consider outflows from rapidly rotating stars in the 
propeller regime. In earlier work we performed multiple sim- 
ulation runs of the propeller stage at a wide variety of pa- 
rameters (R05; U06). Here we take the reference run shown 
in R05 and perform additional analysis. The parameters are 
P* = 1 day, av = 0.3 and ad = 0.2 (see §2.3.2 for the other 
parameter values) . 

Fig.|8]shows snapshots of the matter flow in the propeller 
regime at different times T. One can see that the outflow ap- 
pears at T ^ 50 and continues for a long time (T = 2200 ro- 
tations, or about 6 years in appHcation to protostars). These 
simulations are about 10 times longer than previous simula- 
tions of outflows from a real disk (e.g. Goodson et al. 1997; 
Matt et al. 2002; Kiiker et al. 2003). They are comparable in 
length with simulations of outflows from the disk as a bound- 
ary condition (e.g., Fendt 2009), although here we consider 
outflows form the "real" cold disk to a hot low-density corona. 

Fig. [s] also shows that the outflow has two components. 
One is a conical-shaped wind similar to the conical winds of 
slowly rotating stars discussed earlier. The other component 
is a fast flow of matter interior to the conical winds, which 
we term the axial jet. The disk-magnetosphere interaction 
is strongly non-stationary; the magnetic field lines episodi- 
cally inflate and the disk oscillates. The conical wind com- 



ponent seems to be weakly collimated inside the simulation 
region. The jet component has stronger collimation. The jet 
collimation is stronger in the flow closer to the axis, and 
is enhanced during periods of strong inflation, like at times 
T — 560 and 870. See animations of propeller-driven outflows 
at http://www.astro.cornelledu/^romanova/propellerhtm . 

Fig. [9] shows a typical snapshot from our simulations at 
time T = 1400, with the dimensionless density and velocity 
at sample points (see Table 1 for reference values). One can 
see that the velocities in the conical wind component are sim- 
ilar to those in conical winds around slowly rotating stars. 
Matter launched from the disk has a velocity that is mainly 
azimuthal and approximately Keplerian. It is gradually accel- 
erated to poloidal velocities Vp ^ (0.3 — 0.5) i;^. The flow has 
a high density and carries most of the disk mass into the out- 
flows. The situation is the opposite in the axial jet component: 
the density is 10^ — 10^ times lower, while the poloidal and to- 
tal velocities are much higher. Thus we find a two-component 
outflow: a dense, slow conical wind and a low-density, fast 
axial jet. 

Fig.[To]shows the time-variation of the equatorial density 
and the poloidal magnetic field in the inner part of the sim- 
ulation region. One can see that the density and the poloidal 
magnetic field are strongly enhanced at the inner edge of the 
disk, and the inner disk radius shows large oscillation (see 
also R05, U06). 

Fig. \TT\ l shows the projection of different forces onto a 
closed field line which starts in the disk at r = 4.3 where the 
base of the conical wind (see Fig.|9|. We take only the part of 
the line from the disk to the neutral point where Br = (this 
is the analog of part AB of the line in Fig. [sj) . One can see 
that the forces are large but more or less compensate each 
other. The magnetic force (M) seems to drive matter from 
the disk into the conical wind, though other forces, such as 
the centrifugal (C) and pressure gradient (P) forces are also 
important. It is interesting that conical winds in slowly rotat- 
ing stars and in stars in the propeller regime are similar, but 
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Figure 12. Top panels: matter and energy fluxes atT = 500. Bottom panels: angular distribution of different quantities at radius r = 6 (starting 
from the axis). Panel a: the background and vectors show matter flux Fm, lines are poloidal field lines. The thick dashed line shows the neutral 
hne of the magnetic flux where Br = 0. Panel b: the background shows the energy flux carried by the magnetic field, Fgf, vectors are poloidal 
velocity vectors. Panel c: the background shows the energy flux carried by the matter, Fgf . The solid red line is the line where Vp = va, red 
dashed line corresponds to Vp = i;fm, the solid white line corresponds to Vp = Vsm- Panel d: distribution of density p and matter flux Fm at 
r = 6; Panel e: same but for Fgf and velocities Vp,V(f)', Panel f: same but for Fehi and entropy S. Dashed vertical line shows position of the neutral 
point, Br = 0. 
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Figure 13. The figure is similar to Fig. [12] but different quantities are shown. Panel a: the background shows the distribution of the kinetic beta 
parameter (3i. The solid white line is = 1 line, the red line is (3 = 1 line. Arrows are velocity vectors. Panel h: the background shows angular 
momentum flux carried by the magnetic field, , streamlines show direction of this flux; the thick dashed line corresponds to Br = 0. Panel c: 
the background shows angular momentum flux carried by the matter, FLm- The solid red lines show contours of the flux; Panel d: angular 
distribution of f3 and f3i atr = 6. The dashed vertical line shows the position of the neutral Une Br = 0. Panel e: same but for . Panel f: same 
but for and the viscous flux Flv (dashed line). 
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that the distribution of forces is somewhat different. In coni- 
cal winds the winding of the field lines gives rise to a magnetic 
force in one localized region (above the inner disk) and this 
force dominates. In the propeller regime the disk oscillates 
strongly, and it is important that the magnetosphere presents 
a centrifugal barrier for this matter, and therefore the cen- 
trifugal and pressure gradient forces have a larger role. The 
magnetic force remains important. 

Panel b shows the forces along the coronal field line 
which starts on the surface of the star. We consider the sec- 
ond line from the axis in Fig. [9] which is strongly inflated and 
is a representative line for the description of matter flow into 
the axial jet. One can see that the magnetic force M is much 
larger than the other forces and is the main force accelerating 
matter into the jet. 

Panel c shows velocities along the disk field line (as in 
panel a) . One can see that the azimuthal component dom- 
inates, while the poloidal velocity Vp increases gradually from 
a very small value near the disk up to values comparable with 
Vd). It crosses the slow magnetosonic surface just above the 
disk, and later, the Alfven and the fast magnetosonic surfaces. 

Panel d shows that in the coronal region, the velocities 
are high and the poloidal velocity dominates. Matter crosses 
the slow magnetosonic surface but stays sub-Alfvenic. Both 
the Alfven, va, and the fast magnetosonic, Vfm, velocities are 
about 10 times larger than the flow velocity in the axial jet 
(note the scale at the right-hand side). The flow is in the 
Poynting flux regime found in simulations by Ustyugova et al. 
(2000) and analyzed theoretically by Lovelace et al. (2002). 



4 ANALYSIS OF FLUXES: MATTER, ENERGY, ANGULAR 
MOMENTUM 



4.1 Fluxes in Conical winds 



Fig. \T2\ l shows the matter flux distribution Fm and a neutral 
line of the magnetic field where Br = 0. This line separates 
the field lines starting on the disk from those starting on the 
star. One can see that the conical wind flows along both sets 
of field lines. Panel d shows that the matter flux Fm has a 
sharp peak in its angular distribution (at r = 6), that is, the 
wall of the cone is narrow. The position of the Br = line in 
this panel shows that matter flows along both the stellar and 
the disk field lines. For > 75°, the matter flux is dominated 
by the disk and is negative. We also see that the density p is 
low in the corona. In the disk the density increases to much 
larger values, p = 1 — 10. There is also a low-density gap at 
^ ^ 45° where matter is accelerated to high velocities. 

The energy flux Fe is the sum of the matter component 
Fehi and the field component FEf where 



E X B) 



(5) 



where w is the enthalpy and is the gravitational potential. 

Fig. [T2J) shows that the magnetic field energy flux Fm is 
high near the star, at the base of the conical wind and in the 
region of fast flow. Panel e shows that the energy flux FEf at 
r = 6 has a peak in the region of the conical wind. 

Fig. [T2| : shows that the distribution of the matter en- 
ergy flux, Feui, is similar to the matter flux distribution. The 
panel also shows that in the conical wind, matter crosses 



all the critical surfaces. It ends up flowing with a super- fast- 
magnetosonic velocity. In contrast, in the corona, away from 
the regions of outflow, the flow is sub-slow-magnetosonic. 
Panel / shows that the entropy S is high in the corona and 
low in the conical wind and the disk. 

Fig. [Tsji and d shows the ratio of the gas and magnetic 
pressures, which is the conventional f3 parameter. We also use 
what we term the kinetic parameter /3i, where 
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(6) 



The flow region is magnetically dominated when the (3 or f3i is 
less than unity. The magnetic pressure dominates only in the 
region near the star in the conical wind case. The situation 
is different in the propeller regime where the axial region is 
magnetically dominated (see §4) . 

We calculate the angular momentum flux distribution 
which consists of three components, Fl = Flhi + FLf + Flv, 
where FLm, FLf and Flv are the angular momentum fluxes 
carried by the matter, the magnetic field, and the viscosity: 
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(7) 



Fig. [I3p shows that the magnetic component i^Lf of the flux, 
dominates near the star and in the part of the conical wind 
close to the disk. The streamlines show that angular momen- 
tum flows from the disk onto the star, from the disk into the 
conical wind, and from the star into the corona. Panel c shows 
that the conical winds carry away angular momentum as- 
sociated with matter, Flhi. The magnetic component is also 
high at the base of conical wind. However, at larger distances 
this angular momentum is converted into angular momen- 
tum carried by matter. Comparison of panels e and / shows 
that at r = 6, the angular momentum carried by the mat- 
ter is much larger than that carried by the field. Panel e also 
shows that some angular momentum flows into the disk wind 
(at ^ > 60°). Panel /shows that the angular momentum car- 
ried by viscosity is significant. The disk viscous component is 
much larger than the matter component flowing into the con- 
ical wind, so most of the angular momentum flows outward 
along the disk. 

We also calculated the matter and angular momentum 
fluxes flowing through the surface of the star, and through a 
spherical surface of radius r = 6. 



M = J dS-Fm , Fm = pvp 

L = J dS- (FLm + FLf + Flv) , 



(8) 
(9) 



where dS is the surface area element directed outward. Fig. 
14 (left panels) shows that about 2/3 of the incoming disk 
matter flows to the star, and the rest going into the conical 
wind. Fluxes into the wind oscillate due to magnetic field in- 
flation and reconnection events. The right-hand panels show 
that angular momentum flows inward with the disk matter, 
Ldm. Part of this angular momentum is carried away by the 
conical wind. Angular momentum is carried mainly by mat- 
ter, Lwm >> i^wf. Part of the disk angular momentum flows 
to the star and spins it up (the star rotates slowly in that 
Tcor = 3). The angular momentum carried to the star by 
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Figure 14. Panel a: matter flux onto the star Ms, into the conical wind Mw, and through the disk Md. Panel b: angular momentum flux carried 
by the disk Ld, onto the star carried by matter Lsm and by the magnetic field Lgf • The angular momentum flux of the conical wind carried by 
matter is Lwm and that carried by the field is Lwf • Panels c, d show the time dependence of these fluxes. 



the matter, Lsm, is converted into angular momentum car- 
ried by the field, Lsf, and hence on the surface of the star 
(see also R02). One can see from Fig. 



14 



that all 



Lsf ^sm 

these fluxes are smaller than the flux carried by the disk mat- 
ter. This means that the main part of the angular momentum 
of the disk flows outward to larger distances due to viscosity. 



4.2 Fluxes in the propeller regime 

We analyze fluxes in the propeller regime in a similar manner 
to that for the conical winds. Fig.[T5|i shows the distribution 
of the poloidal matter flux Fm in the background with the 
vectors Fm on top. The neutral line (dashed white line) sep- 
arates the field lines which start in the disk from those which 
start on the star. 

Panel d shows that most of the matter flows along field 
lines threading the disk, although some matter flows along 
the stellar field lines. For ^ > 80° (in the disk region) the mat- 
ter flux becomes much larger and negative (we exclude this 
part of the plot to show the conical wind part more clearly) . 
The plot of the density shows that the density is very low on 
the axis but gradually increases toward the region of the con- 
ical wind and continues to grow towards the disk. 

Panel b shows the distribution of the magnetic energy 
flux FEf. One can see that a strong flux of magnetic energy 
(Poynting flux) flows into the corona. This is the region where 
matter is accelerated to high velocities (see velocity vectors) . 
Panel e shows that at r = 10 the magnetic energy flux is very 
large, and is distributed over a range of angles with a maxi- 
mum at ^ 25° (not on the axis). The plot also shows that 
the poloidal velocity is slightly larger than the azimuthal ve- 
locity and both velocities are high, up to = 2 (which is 400 
km/s for protostars). The jet component is smaller in the case 
of slowly rotating stars. 

Panel c shows the energy flux associated with the matter 



flow. One can see that matter in the conical wind crosses the 
slow magnetosonic surface v — Vsm just above the disk, and 
soon crosses the Alfven surface v = va, and the fast magne- 
tosonic surface v = Vim- Panel /shows that the matter energy 
flux distribution has a sharp peak in the region of the conical 
wind and that the entropy S is high in the corona but drops 
towards the disk. 

Comparison of panels e and/ shows that the maximum of 
the energy flux carried by the magnetic field into the jet, FEf, 
is about 3 times larger than that carried by the matter into 
conical winds, FEm- In addition, the integrated flux carried 
by the magnetic field is a few times larger. Therefore, the cu- 
mulative energy flux carried by the magnetic field into the jet 
(the Poynting flux) is about 10 times larger than that carried 
by matter. This means that the jet component is 10 times more 
powerful. Part of this energy is converted into kinetic energy 
of the fast component inside the simulation region. However, 
most of the magnetic energy may be transferred to particles 
or converted into radiation at larger distances from the star. 

Next, we analyze the angular momentum flow. Fig. [T6| i 
shows the /3-parameter in the background and the (3 = 1 and 
/3i = 1 surfaces in the foreground (see equation 6). One can 
see that the magnetic energy dominates in the whole axial 
coronal region interior to the conical wind (/3i < 1 and [3 < 
1). If one uses only the standard criterion /3 = 1, then one 
can see that the region above the disk is also magnetically 
dominated (/3 < 1). Panel d shows the angular distribution of 
(3 and /3i at r = 10. 

Panel b shows the distribution of the angular momentum 
flux carried by the magnetic field, FLf , and the streamlines as- 
sociated with this flux. One can see that a significant amount 
of angular momentum flux flows from the star into the corona 
along the stellar field lines. Panel e shows that the angular 
momentum flows out along the set of field lines between the 
axis and the neutral field line with a maximum right above 
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the conical wind component. Some of them thread the low- 
density corona, while others thread the upper part of the con- 
ical wind above the neutral line. These panels also show that 
a significant amount of angular momentum flows from the 
inner part of the disk. See U06 for detailed analysis of the 
different components of the angular momentum flow. 

Panel c shows that the angular momentum flux carried 
by the matter is also large and is carried by the conical wind. 
Panel / shows that most of this angular momentum flows 
along the disk field lines while some angular momentum 
flows along the stellar field lines. 

Fig. [iTp shows the matter fluxes onto the star Ms, and 
into the outflows, M^, integrated over a surface with radius 
r = 10 (any flow with > is taken into account). One 
can see that the matter flux into the wind is much larger 
than that onto the star, Mu, » Ms, that is, almost all disk 
matter is ejected from the system into the outflows. Here 
we should note that we consider the "strong propeller" case, 
Vcor « Tm- If thc Star rotates slower, then the fraction of 
the matter flux going into the wind decreases, and a larger 
portion of the matter may accrete onto the star (see U06 for 
dependences of matter fluxes on Q*, B^, av and a^.) Both 
fluxes are strongly variable and show episodic enhancement 
of accretion and outflows. Simulations show that an interval 
between the strongest outbursts increases when diffusivity co- 
efficient ad decreases (R05, U06). 

Panel h shows the integrated angular momentum fluxes 
through the same r = 10 surface. Here we calculate sepa- 
rately the angular momentum fluxes carried by the field and 
by the matter. One can see that the star spins down due to the 
angular momentum carried by the magnetic field, Lsf, while 
the angular momentum carried by the matter flow, Lsm, is 
negligibly small. The angular momentum outflow from the 
star, Lsf, almost coincides with the angular momentum car- 
ried by magnetic field lines into the magnetically-dominated 
jet, Lwf. This indicates that angular momentum flows from 
the star into the magnetically dominated axial jet. Thus a star 
in the propeller regime is expected to spin down rapidly due 
to angular momentum flow into the magnetically-dominated 
axial jet. Analysis of U06 shows that this jet angular momen- 
tum is approximately equally split between the flux carried 
into the corona along open field lines, and the flux which 
flows along partially inflated field lines which close inside the 
simulation region and are connected with both the star and 
the disk. Panel c shows that the angular momentum carried 
by matter into the conical winds, Lwm is approximately equal 
to that carried by the field to the corona. The bottom panels 
show the same plots at higher time resolution. Therefore, the 
star-disk system loses its angular momentum through both 
the wind and jet components, via the inner disk and star re- 
spectively. So, there is no problem with excess angular mo- 
mentum in the star-disk-system; it flows into the jet/wind. 



5 OTHER PROPERTIES OF OUTFLOWS 

5.1 Inflation of Field Lines and Disk Oscillations 

The field lines connecting the disk and the star have the 
tendency to inflate (e.g., Lovelace, Romanova & Bisnovatyi- 
Kogan 1995). Quasi-periodic reconstruction of the magneto- 
sphere due to inflation and reconnection has been discussed 



theoretically (Aly & Kuijpers 1990; Uzdensky Litwin & Konigl 
2003) and has been observed in a number of axisymmetrtic 
simulations (Hirose et al. 1997; Goodson et al. 1997, 1999; 
Matt et al. 2002; Romanova et al. 2002 - hereafter R02; von 
Rekowski & Brandenburg 2004). Goodson & Winglee (1999) 
discuss the physics of inflation cycles. They have shown that 
each cycle of inflation consists of a period of matter accumula- 
tion near the magnetosphere, diffusion of this matter through 
the magnetospheric field, inflation of the corresponding field 
lines, accretion of some matter onto the star, and outflow of 
some matter as winds, with subsequent expansion of the mag- 
netosphere. There simulations show 5 — 6 cycles of inflation 
and reconnection. 

Our simulations show 30 — 50 cycles of inflation and re- 
connection in the propeller regime. We chose one outburst 
from our simulations and plotted the density and a fixed set 
of magnetic field lines at different times. Fig. [Ts] shows that at 
T — 890, the magnetosphere is relatively expanded, although 
some matter accretes around the expanded field lines (see 
also Romanova et al. 2004a). At T = 900, the disk matter 
comes closer to the star and some field lines inflate or par- 
tially inflate, thus blocking accretion. At T = 910 even more 
field lines inflate and accretion is blocked. However, outflow 
is permitted at both of these moments of time. At T = 916, 
the internal field lines reconnect, permitting accretion onto 
the star. At T = 921, the magnetosphere expands and accre- 
tion onto the star is again prevented. Later, at T = 927, the 
field lines reconnect and some matter accretes along a longer 
path - around the expanded magnetosphere. This picture is 
similar to that described by Goodson & Winglee (1999). 

Fig. [TTji shows that the time interval between the 
strongest outbursts in the propeller regime is AT 50 — 70. 
In appHcation to protostars and CTTSs (Po = 1-04 days) this 
time corresponds to (At)outb = 52 - 73 days. In some young 
stars, like CTTS HH30 (XZ Tau), for example, the outbursts 
into the jet occur at intervals of a few months, which hints 
that episodic inflation of field lines may be responsible for 
some outbursts. During the outbursts, the matter flux into 
outflows increases several times and the velocities also in- 
crease. This may lead to the formation of new blobs or to the 
generation of shock waves in the outflow. This mechanism 
may be relevant for formation of blobs or shocks in protostars 
and rapidly rotating CTTSs. In slowly rotating stars, the time- 
interval between outbursts is smaller, (At)outb ~ 5 days (see 
Fig. [T4I ), so the outbursts have a smaller amplitude but are 
more frequent. The interval depends on the diffusivity in the 
disk, ttd. At very small diffusivity the time-interval between 
outbursts may be much larger. 

Diffusivity is important for reconnection processes in the 
corona. We have diffusivity only in the disk. We choose a 
certain density level pd — 0.3 below which the diffusivity 
is absent, so that high-density regions, p > 0.3, which cor- 
respond to the disk and the funnel streams, have diffusiv- 
ity, and low-density regions do not. In the corona and the 
conical outflows, the diffusivity has only a numerical origin 
and is small. Namely, we observe in simulations that in con- 
ical winds, the two layers of plasma with an oppositely di- 
rected magnetic field reconnect only slowly. Similar behavior 
has been observed in ideal MHD simulations by Fendt & El- 
stner (2000) . An anomalous (high) diffusivity was added by 
Hayashi et al. (1996) to a part of the simulation region to 
enhance the reconnection process in the inflating plasmoids. 
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Diffusivity had been added into the whole simulation region 
by Fendt & Cemeljic (2002). They observed that at higher dif- 
fusivity the level of collimation by the magnetic field and the 
Lorentz force decrease, while the centrifugal force increases. 
We performed exploratory simulations with non-zero diffusiv- 
ity in the corona. We added to the corona the same diffusivity 
as in the disk with = 0.1 — 0.2, which operates at different 
density levels, p > pd, where pd = 0.1,0.03,0.01 and in the 
whole simulation region (formally, p^z = 0) . We observed that 
in case of conical winds (slowly rotating stars) the diffusivity 
in the corona does not change the result. However, in case of 
propeller-driven winds, we observed that propeller becomes 
weaker. We believe that the difference is in the fact that in 
the case of conical winds, the wind and the neutral line of 
the inflated magnetic field have approximately the same po- 
sition in space, leading to slower reconnection. In the other 
case, in the propeller regime the inner disk and the region of 
outflows strongly oscillates, and so the position of the neu- 
tral line varies, and hence the reconnection is forced (that is, 
the plasma layers with oppositely directed fields are pushed 
towards each other by an external force) . 

5.2 Matter loading onto stellar field lines, and possible 
role of stellar wind 

Here we discuss how the disk matter gets loaded onto the 
stellar field lines and then flows into the jet in the propeller 
regime (where the jet is strong). It is important to have dif- 
fusivity in the disk, so that the matter of the disk threads the 
field lines of the star and flows onto the star in funnel streams. 
When a sufficient amount of matter is accumulated in the in- 
ner disk, the field lines connecting the star and the disk in- 
flate. During and after inflation, part of the disk matter ends 
up on the field lines connecting the disk with corona (usu- 
ally most of the matter flows along these field lines) . Another, 
smaller part of the matter ends on the field lines connecting 
the star with the corona. For example. Fig. 15a demonstrates 
the result of such inflation, where the neutral line dividing 
the stellar and disk lines is in the middle of the conical wind 
component. On the other hand, when matter flows in a funnel 
stream, most of it accretes onto the star. However, part of it 
is stripped away by the magnetic and centrifugal forces and 
flows into the jet along the stellar field lines. Fig. [19] shows 
that there is a dividing line running through the upper part of 
the funnel stream, separating the regions from which matter 
flows onto the star (most of it) from those from which it flows 
into the jet along the stellar field lines (a small fraction). In 
the funnel region the density is usually high enough so that 
the diffusivity which works in our disk also works in the fun- 
nel stream. This diffusivity helps launch matter from the fun- 
nel stream field lines onto the coronal, jet field lines. Both 
processes are consistent with the strong decrease in coronal 
matter density along the axis. This region is "matter-starved". 

Our simulations do not take into account possible stellar 
winds. Even a weak wind from the star may have a significant 
influence on the axial region of the jet in the propeller regime 
and the "matter-starved" jet region in slowly rotating stars. 
The existence of powerful stellar winds was suggested by Matt 
& Pudritz (2005, 2008) in order to explain the loss of angular 
momentum by young stars. The spectra of many CTTSs (e.g. 
Edwards et al. 2003; Dupree et al. 2005) require up to 10% 
of the disk mass flowing out as winds, in order to explain 



different spectral lines (Edwards 2009). No such winds are 
observed in diskless, weak-line T Tauri stars. Hence, the winds 
must be accretion-driven (e.g., Edwards et al. 2006; Kwan et 
al. 2007). The physics of these accretion-driven stellar winds 
is not understood yet. In the standard approach it is suggested 
that matter falling onto the surface of the star through the 
funnel stream forms a shock near the surface and is heated 
by this shock. However, it cools rapidly in the radiative zone 
behind the shock wave, and no reverse flow into the wind is 
expected (e.g., Lamzin 1998; Koldoba et al. 2008). In another 
investigation, however, Alfven waves and other processes at 
the stellar surface help accelerate up to 1% of the accreting 
disk matter into the wind (e.g. Cranmer 2008). We did not 
incorporate stellar winds into the present simulations. Weak 
winds may help supply matter to the magnetically-accelerated 
axial jets and the "matter-starved" region of fast flow in slowly 
rotating stars. On the other hand, if the wind is very strong, 
say Msw O.lMd, it will probably be matter-dominated at 
moderate distances from the star (say > 10 stellar radii) and 
will have a decollimating effect on the outflows (Fendt 2009) . 
In summary, a weak stellar wind will contribute matter to the 
jet component. 

5.3 Collimation of outflows 

Collimation of conical winds. We observe conical winds in 
both slowly and rapidly rotating stars. In both cases, matter in 
the conical winds passes through the Alfven surface, beyond 
which the flow becomes matter-dominated. We note that in 
slowly rotating stars, the distribution of the poloidal current 
Jp (see Fig.js]:) is such that the corresponding magnetic force 
has a component towards the axis. This may explain why 
conical winds show some collimation (see Fig. [2| . The con- 
ical wind component of the propeller-driven outflows shows 
stronger collimation during periods of inflation and outbursts 
(see Fig. [8| . However, this collimation may not be sufficient 
to explain well-collimated jets. 

Conical winds may be further collimated at larger distances 
from the star either by the pressure of the external medium 
(Lovelace et al. 1991; Frank & Mellema 1996), or by disk 
winds (Konigl & Pudritz 2000; Ferreira et al. 2006; Matsakos 
et al. 2008; Fendt 2009). In addition. Matt, Winglee, & Bohm 
(2003) have shown that a weak axial magnetic field {B « 
0.1 G) associated with the disk, may collimate the winds at a 
distance of a few AU. 
Collimation of the jet. In the propeller regime, the jet com- 
ponent is self-collimated by the magnetic hoop-stress. The 
level of collimation increases towards the axis. The poloidal 
velocity in the jet also increases towards the axis, and varies 
between Vp ^ 2 near the axis and Vp ^ 0.2 near the coni- 
cal wind. That is why we choose a few typical velocity levels 
Vp = Vc, with Vc = 0.5, 1 and 1.5, and plot lines of equal 
velocity (Fig. [20] ) . In application to protostars, the fast com- 
ponent of the jet, Vp > 200 km/s, carries ^ 2% of the mass 
and ^ 22%/0.6 ^ 37% of the angular momentum flux out 
of the star. At the lower velocity limit, Vp > 100 km/s, these 
numbers are 10% and 35%/0.6 ^ 60%. 

It is of interest to know the dependence of the mass out- 
flow rate on the poloidal velocity Vp. We calculate the matter 
flux M{vp > Vc) through the external boundary r = Rout at 
poloidal velocities above a certain value Vc for different val- 
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Figure 15. Same as in Fig.[T2]but for the propeller regime at time T = 1400. 
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Figure 16. Same as in Fig.fTslbut for the propeller regime at time T = 1400. 



ues of Vc. Panel a shows that the jet component with Vp > 0.5 
carries about 10% of the total outflowing mass Mu,, while the 
very fast components, Vp > 1.5 and Vp > 1 carry only about 
1% and 2% correspondingly. So, about 10% of the total mass 
flows into the collimated jet. The fractions of matter flowing 
into different parts of the jet and into the conical wind are 
shown in the left panel. 

The jet carries angular momentum out of the star along 
different field lines corresponding to different Vp. It is of inter- 
est to know which part of the jet carries most of the angular 
momentum. We calculate the angular momentum flux carried 



by the magnetic field Lf{vp > Vc) (this component dominates 
in the jet) through the external boundary and normalize it to 
the total magnetic flux through this boundary Lwf. Panel b 
shows that the high-velocity part of the jet, Vp > 1.5, carries 
about 13% of the total angular momentum flux, while the en- 
tire inner part of the jet in the velocity intervals Vp > 1 and 
Vp > 0.5 carry 22% and 35% of the flux correspondingly. An- 
other fraction (12%) flows out along stellar field lines thread- 
ing the conical wind component and the low-velocity area 
above it. All this flux is responsible for spinning down the 
star (which is about 60% of the total flux). The rest of the 
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Figure 17. Panel a: matter fluxes onto the surface of the star Ms, and into the outflows (both wind and jet) through the surface r = 10 versus 
time T. Panel b: angular momentum fluxes carried to (or out from) the star by the matter, Lgm (gray line), and by the magnetic field, Lgf (red 
hne). The green line shows the angular momentum flux carried by the magnetic field through the surface r = 10. Panel c: angular momentum flux 
carried by the matter component Lwm through the surface r = 10. The bottom panels show the same fluxes but during a part of the simulation 
time. 
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Figure 18. Several snapshots from simulations show the mechanism of enhanced accretion and outflows observed in the light-curves of propeller- 
driven accretion and outflows. The background shows the matter flux, lines are the sample magnetic field lines, the same set of which is shown 
in all plots. 



flux (40%) flow^s along the disk field lines threading the con- 
ical winds and the disk. We conclude that the jet component 
above the conical wind carries a relatively small mass but has 
a significant contribution to the angular momentum outflow 
from the star. Note that only about half the star's angular mo- 
mentum flows into the jet. The other half is associated with 
star-disk interaction through the field lines which are closed 
inside the simulation region and were not taken into account 
in this analysis (see U06 for details). In application to proto- 
stars, the fast component of the jet, Vp > 200 km/s, carries 
- 2% of the mass and - 22%/0.6 ^ 37% of the angular 
momentum flux out of the star. At the lower velocity limit, 
Vp > 100 km/s, these numbers are 10% and 35%/0.6 ^ 60%. 




Figure 19. White lines with arrows show streamlines of matter flow 
in the propeller regime at T = 1500. Red lines show sample magnetic 
field lines. Only the inner part of the simulation region is shown. 



Launching of Conical Winds 1 7 



22% I 



f V p cf V. 




1000 



Time 



2000 



1000 Time 2000 



Figure 20. Analysis of matter and angular momentum distribution between the jet and wind components in the propeller regime at T = 1400. 
Left panel: The yellow lines show surfaces of constant poloidal velocity, i;c = 0.5,l,1.5(in application to protostars they are approximately 100, 
200 and 300 km/s). Numbers on top of the plot show fractions of total matter flux Mw into different sectors of the region, and the same for 
angular momentum carried by the field, L^f- Panel b: Matter flux into the jet through an external boundary, r = Rout, with velocities Vp > Vc, 
M{vp > Vc) versus total matter flux Mw Panel c: Same as in panel b, but for the angular momentum flux carried by the magnetic field. 



6 3D SIMULATIONS OF CONICAL WINDS 

We did exploratory simulations of conical winds in global 
3D simulations. We chose a case where the dipole magnetic 
field of the star is misaligned with the rotation axis (of the 
star and disk) by an angle 6 = 30°. One question is what 
the direction of the conical wind is in the case of an in- 
clined dipole. We used the Godunov-type 3D MHD "cubed 
sphere" code developed by Koldoba et al. (2002) . In the past 
we have used this code to study magnetospheric accretion 
close to the star (Romanova et al. 2003, 2004b). Compared 
with that work we decreased the density in the corona by a 
factor of 10 to pc = 0.001 and created conditions suitable 
for bunching of the field lines. We used a grid resolution of 
AT^ X AT^ = 120 X 51^ in each of 6 blocks of the sphere. We 
took the density in the disk to be = 2 which is 5 times lower 
than in the axisymmetric case shown above. At the same time, 
we chose a smaller magnetic moment for the star, /x = 2 com- 
pared to /i = 10 in the axisymmetric case (to reduce the com- 
puting time) . We start the disk flow not from large distances 
but from r = 5, to limit the computing time. The bunching 
of field lines is achieved by having a sufficiently high viscos- 
ity, av — 0.3. We do not have diffusivity in the 3D code, but 
at the grid resolution we use, the estimated numerical diffu- 
sivity at the disk-magnetosphere boundary is at the level of 
ad ^ 0.01 — 0.02, and hence the main condition for conical 
wind formation, Prm ^ 1, is satisfied (see also Appendix C). 

Simulations show that the accreting matter bunches up 
field lines and some matter flows out as a conical wind. Fig. 
[2T] shows that the wind is geometrically symmetric about the 
rotation axis. However, the density distribution in the wind 
shows a spiral structure which rotates with the angular ve- 
locity of the star, Q*, and represents a one-armed spiral wave 
from each side of the outflow. 

Note that for high av (0.3 in this case), the disk- 
magnetosphere boundary may exhibit the magnetic inter- 
change instability (Romanova, Kulkarni & Lovelace 2008; 
Kulkarni & Romanova 2008) . In these simulations we do ob- 
serve some accretion due to this instability in addition to the 
main funnel stream accretion that dominates at high mis- 
ahgnment angles, 6 > 30° (Kulkarni & Romanova 2009). 



However, the conical wind originates at larger radii compared 
with the inner disk radius where accretion through instability 
dominates. We believe that both processes can "peacefully" 
co-exist for 6 > 30°. However, in other situations the coni- 
cal wind may be influenced by the interchange instability. For 
example, we did not try to investigate outflows at small B 
where accretion through instability often dominates. Accre- 
tion through instability opens up a new path for penetration 
of matter through the magnetosphere, and thus may possi- 
bly decrease the bunching of field lines and consequently the 
strength of conical winds. This interrelation between instabil- 
ities and conical winds needs to be investigated in future 3D 
simulations. Longer simulations should be performed, and ac- 
cretion to rapidly rotating stars should also be examined. 



7 COMPARISON WITH THE X-WIND MODEL 

Winds from the disk-magnetosphere boundary have been pro- 
posed earlier by Shu and collaborators and referred to as X- 
winds (e.g., Shu et al. 1994). In this model, X-winds originate 
from a small region near the corotation radius rcor, while the 
disk truncation radius n (or, the magnetospheric radius rm) 
is only sHghtly smaller than rcor (rm ~ 0.7 rcor, Shu et al. 
1994). It is suggested that excess angular momentum flows 
from the star to the disk and from there into the X-winds. The 
model aims to explain the slow rotation of the star and the 
formation of jets. In the simulations discussed here we have 
obtained outflows from both slowly and rapidly rotating stars. 
Both have conical wind components which are reminiscent of 
X-winds. What, then, is the difference between X-winds, con- 
ical winds and propeller-driven winds? 

In some respects conical/propeller winds are similar to 
X-winds: (1) They both require bunching of the poloidal field 
lines and show outflows from the inner disk; (2) They both 
have high rotation and show gradual poloidal acceleration 
(e.g., Najita & Shu 1994). 

The differences are the following: (1) The coni- 
cal/propeller outflows have two components: a slow high- 
density conical wind (which can be considered as an ana- 
logue of the X-wind), and a fast low-density jet. No jet com- 
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ponent is discussed in the X-wind model. (2) Conical winds 
form around stars with any rotation rate including very slowly 
rotating stars. They do not require fine tuning of the coro- 
tation and truncation radii. For example, bunching of field 
lines is often expected during periods of enhanced or unsta- 
ble accretion when the disk comes closer to the surface of the 
star and rm « Vcor- Under this condition conical winds will 
form. In contrast, X- winds require rm ~ Vcor- (3) The base 
of the conical wind component in both slowly and rapidly ro- 
tating stars is associated with the region where the field lines 
are bunched up, and not with the corotation radius. (4) X- 
winds are driven by the centrifugal force (Blandford & Payne 
1982), and as a result matter flows over a wide range of direc- 
tions below the "dead zone" (Shu et al. 1994; Ostriker & Shu 
1995). In conical winds the matter is driven by the magnetic 
force (Lovelace et al. 1991) which acts such that the matter 
flows into a thin shell with a cone angle ^ 30°. The same 
force acts to partially collimate the flow. (5) In the X-wind 
model it is suggested that angular momentum flows from the 
star to the disk in spite of the fact that the truncation radius 
of the disk is located at rm ~ 0.7rcor and the disk rotates 
faster than the star (Shu et al. 1994). Simulations show that 
if the funnel stream starts at Vm < rcor, then angular mo- 
mentum flows from the disk to the star along magnetic field 
lines of the funnel stream which form a leading spiral, and 
the star spins up (R02, Romanova et al. 2003; Bessolaz et al. 
2008). The star may transfer its angular momentum to the 
disk if rm > rcor, like in the propeller case considered above. 
(6) The X-wind regime is somewhat similar to the propeller 
regime, where the star transfers part of its angular momen- 
tum to the disk, and this excess angular momentum may flow 
into the conical component of the wind. However, in the pro- 
peller regime, angular momentum also flows from the star 
into the jet. (7) Conical and propeller-driven winds are non- 
stationary: the magnetic field constantly inflates and recon- 
nects. X- winds, on the other hand, are steady. This difference, 
however, is not significant, and models can be compared us- 
ing time-averaged characteristics. 



8 APPLICATION TO DIFFERENT STARS 
8.1 Application to Young Stars 

Our simulation results can be applied to different types of 
young stars, including low-mass protostars (class I YSOs) 
which often show powerful outflows, CTTSs (class II YSOs) 
which show less powerful outflows, EXors which show peri- 
ods of strongly enhanced accretion and outflows, and young 
brown dwarfs. 



8.1.1 Low-mass protostars (class I YSOs) 

Class I protostars are young stars which are usually embed- 
ded inside a cloud of gas and dust. IR observations show that 
protostars are surrounded by cold massive disks and that the 
accretion rate is usually an order of magnitude larger than in 
CTTSs, that is, M - (10"^ - lO"'^)M0/yr (e.g., Nisini et al. 
2005). The outflows are also more powerful than in CTTSs. 
The stars are fully convective, and so rapid generation of a 
magnetic field that may even be larger than in CTTSs is ex- 
pected. 
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Figure 21. Projections show formation of conical winds obtained in 
3D MHD simulations. Left panels show the distribution of matter flux 
(background) and sample magnetic field lines in the XZ plane at dif- 
ferent moments of time. Right panels show corresponding YZ slices. 
Arrows show the direction of the magnetic moment of the star /i and 
angular velocity of rotation . 



We consider a protostar of mass Mo = M = O.SMq, ra- 
dius = 2i?0, and surface magnetic field = 3 x 10^ G. 
The dimensionless radius of the star (the inner boundary) is 
0.5, and the unit radius Ro = 2i?* = 2.8 x 10^^ cm. The veloc- 
ity scale is vo = {GM/RqY^'^ = 195 km/s, the time-scale is 
to = i?o/'^^o = 0.16 days, and the period of rotation at i?o = 1 
is Po = 1.04 days. We take a rapidly rotating star with period 
P* = 1.04 days (the corotation radius of Vcor = 1). The other 
reference variables are shown in Table 1. For dimensionless 
temperatures in the disk and corona ofT^ = 5 x lO"'^ and 
Tc — 0.5, we obtain corresponding initial dimensional tem- 
peratures: Td = 2290 K and Tc = 2.3 x 10^ K. Fig. 
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shows 

the distribution of density and velocity around the protostar. 
The age of protostars is 10^ — 10^ years, and therefore they 
may rotate more rapidly than CTTSs and it is likely that some 
of them are in the propeller regime. If the propeller is strong 
enough (like in our simulations, where the period P* 1 day 
and av — 0.3) then most of the disk matter will be ejected as 
slow conical winds with velocity ~ 50 km/s, which may be 
higher if the disk is closer to the star. Most of the energy, how- 
ever, flows into the magnetically-dominated axial jet, where a 
small fraction (about 10%) of the disk matter is accelerated up 
to Vp 100 — 400 km/s inside the simulation region. A huge 
amount of angular momentum flows out of the star through 
the same jet, and conical winds carry a comparable amount 
of angular momentum as well. This may solve the angular 
momentum problem of the system. So, at this stage the out- 
flows are powered by two things: the stellar rotational energy 
and the inner disk winds (the conical winds) . Fig. \\7\ ) shows 
that the outflow is strongly non-stationary with strong matter 
ejection into jets/winds every 2 — 3 months. Ejection is accom- 
panied by larger than average matter flux and velocities, and 
hence formation of new blobs or shock waves is expected. 

A protostar in the propeller regime loses its angular mo- 
mentum to an axial jet. From the right-hand panels of Fig. 
[Tt) we obtain the dimensionless value of the angular mo- 
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mentum loss: Lsw ~ 3, which corresponds to a dimensional 

value of Lsw = LswLo ~ 9.3 x 10^'^ gcm^/s^. The star's 
angular velocity is = 27r/P* !^ 7 x 10~^s~^, its an- 
gular momentum is J = kMr'^Q^ - 2.2 x lO^^k gcm^/s, 
where k < 1. Taking k = 0.4, the spin-down time-scale is 
T = J/ Lsw ~ 3 X 10^ years. Note that this time-scale is cal- 
culated for = 3 X lO^G. The time-scale decreases with the 
magnetic field of the star as ^ B~^-^ (see U06) and will be 
T 3 X 10^ years for = lO^G. If the magnetic field is 
weaker, then the protostar will continue to spin rapidly even 
in the CTTSs stage. U06 present the dependence of the spin- 
down time-scale on the magnetic field, the spin of the star 
and other parameters. 



8. 1 .2 Classical T Tauri Stars (class II YSOs) 

CTTSs and their jets have been extensively studied in recent 
years. High-resolution observations of CTTSs show that the 
outflows often have an "onion-skin" structure, with better- 
collimated, higher-velocity outflows in the axial region, and 
less-collimated, lower-velocity outflows at a larger distance 
from the axis (Bacciotti et al. 2000). In other observations, 
high angular resolution [Fell] A 1.644/im emission line maps 
taken along the jets of DC Tau, HL Tau and RW Aurigae reveal 
two components: a high-velocity well-collimated extended 
component with velocity v ^ 200 — 400 km/s, and a low- 
velocity ^; '-^ 100 km/s, uncollimated component closer to the 
star (Pyo et al. 2003, 2006). High-resolution observations of 
molecular hydrogen in HL Tau have shown that at small dis- 
tances from the star, the flow shows a conical structure with 
outflow velocity - 50 - 80 km/s (Takami et al. 2007). In XZ 
Tau, two-component outflows are observed: one component 
is a powerful but low-velocity conical wind with an opening 
angle of about 1 radian, and the other is a fast well-collimated 
axial jet (e.g., Krist et al. 2008). The origin of these outflows is 
not known, but we can suggest that at least the lower-velocity 
component may be explained by the conical winds suggesting 
that the condition for bunching, Prm > 1, is satisfied. If a 
CTTS rotates rapidly (in the propeller regime) then the jet 
component may originate from the propeller effect. 

Spectral observations of the He I (10830A) line show 
clear evidence of two-component outflows (Edwards et al. 
2003, 2006; Kwan et al. 2007). Observations show (see Fig. 
[23] ) that at smaller accretion rates only the relatively low- 
velocity component, ^; < 100 km/s, appears. At higher ac- 
cretion rates there is evidence of a very fast component with 
V ~ 200-400 km/s which requires an outflow rate of up to 
Mw ^ O.lMd (Edwards et al. 2006; Edwards 2009). We sug- 
gest that the low- velocity component may be a conical wind. 
However, it is not clear what can explain the high-velocity 
component. Even if a star is in the propeller regime, then at 
the required velocities v > 200 km/s, only (2-3)% of the disk 
matter flows into the jet. Possibly, additional matter influx 
from the wind of the star may enhance the matter flux into 
the jet component. In another example, observations of the 
Hp spectral line in RW Aurigae, and comparison of possible 
outflow geometries led to the conclusion that a thin cone- 
shaped wind with a half-opening angle of 30° — 40° gives the 
best fit to the observations (from Alencar et al. 2005; see Fig. 
HI). 

The strongest outbursts supplying CTTS jets are usually 



episodic or quasi-periodic (e.g., Ray et al. 2007). For exam- 
ple, blobs are ejected every few months in HH30 (XZ Tau), 
and every 5 years in DC Tau (Pyo et al. 2003) . Both of these 
may be connected with episodes of enhanced accretion and 
formation of conical winds. The velocity and density in the 
outflow are larger during periods of enhanced accretion, be- 
cause the disk comes closer. If the CTTS is in a binary system, 
then the accretion rate may be episodically enhanced due to 
interaction with the secondary star, and this may explain the 
longer intervals of a few years) between outbursts observed 
in other CTTSs. Events of fast, implosive accretion are also 
possible due to thermal or global magnetic instabilities (e.g., 
Lovelace et al. 1994). Alternatively, a period of a few months 
may be connected with long-term episodes of oscillations of 
the magnetosphere. In the propeller regime the time-interval 
between oscillations is 1-2 months even for mild parameters. 
Bouvier et al. (2007) have shown that magnetospheric ex- 
pansion in the CTTS AA Tau may occur with a period of a few 
weeks. Multi-year observations of variability in CTTSs show 
that they are strongly variable on different time-scales (e.g., 
Herbst et al. 2004; Grankin et al. 2007) which is probably 
connected with periods of enhanced accretion. 

For CTTSs we suggest the same parameters as for pro- 
tostars but take a weaker magnetic field, B = 10^ G, so 
that for the same dimensionless runs we obtain lower ac- 
cretion rates (see Table 1). Taking from Fig. [14J: the dimen- 
sionless values of the matter flux onto the star, Ms ~ 3.8, 
and into the conical winds, Mw ~ 1.3, and taking the value 
of Mo Table 1, we obtain an accretion rate onto the star of 
Ms = MsMo ^ 7.6 X 10"^ Mo/yr and into the wind of 
Mu, ~ 2.6 X 10~^ Mo/yr. For a corotation radius of rcor = 3, 
the period of the CTTS is P* = 5.6 days. In typical simulation 
run the truncation radius rm ~ 1.2 is much smaller than the 
corotation radius rcor = 3. This situation corresponds to the 
case of enhanced accretion when the star spins up, and which 
corresponds to ejection of conical winds. 

Many CTTSs are expected to be in the rotational equilib- 
rium state, when rm ~ rcor - Without the bunching condition, 
and at small viscosity and diffusivity parameters, no signifi- 
cant outflows had been observed in simulations (R02; Long 
et al. 2005). On the other hand, if the bunching condition is 
satisfied and/or the accretion rate is enhanced, then conical 
winds are expected. It is possible that the jet component is 
also powerful enough in this state so as to produce the fast 
jet component that is observed. Additional simulations are 
needed for better understanding of outflows in this important 
state. 



8.1.3 Periods of enhanced accretion and outflows in EXors 

EXors represent an interesting stage of evolution of young 
stars where the accretion rate is strongly enhanced and pow- 
erful outflows are observed (e.g., Coffey, Downes & Ray 2004; 
Lorenzetti et al. 2006; Brittain et al. 2007). Brittain et al. 
(2007) reported on the outflow of warm gas from the inner 
disk around EXor VI 647, observed in the blue absorption of 
the CO line during the decline of the EXor activity. They con- 
cluded that this outflow is a continuation of activity associ- 
ated with early enhanced accretion and bunching of the mag- 
netic field lines (see Fig. [25] ) . The EXor stage may correspond 
to the initial stage of our simulations, during which a signif- 
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Figure 22. A dimensional example of matter flow in the protostar regime shown in Fig. [9] Time T = 1400 corresponds to 3.8 years. Labels show 
the particles density n in units of and the poloidal velocity Vp. Azimuthal velocity is a few times larger that the poloidal velocity in the 

beginning of the flow, but decreases at larger distances (see Fig. |9} . 



icant amount of matter comes into the region. Or, it is more 
probable that initially there is weak outflow at the level of that 
in CTTSs, but later the accretion rate increases by a few orders 
of magnitude, leading to a powerful outburst which produces 
conical winds. For conversion into dimensional values, we 
suggest that the disk comes close to the stellar surface, which 
is at r = 1 (as opposed to 0.5 in the previous examples), and 
the disk stops much closer to the star (rm = 1.2i?*). Then all 
velocities are higher by a factor of V2 ^ 1.4, densities by a 
factor of 32, and matter fluxes by a factor of 11 than in the 
main example relevant to CTTSs. 



8.1.4 Outflows from Brown Dwarfs 

Recently outflows were discovered from a few brown dwarfs 
(BDs) (e.g., Mohanty, Jayawardhana & Basri 2005; Whelan, 
Ray & Bacciotti 2009). Clear signs of CTTS-like magneto- 
spheric accretion (broad spectral lines with full-widths of 
V > 200 km/s) were reported earlier for a number of young 
BDs (e.g., Natta et al. 2004). BDs are fully or partially con- 
vective and the generation of a strong magnetic field is ex- 
pected (Chabrier et al. 2007). Magnetic fields of the order of 
0.1 — 3 kG may explain the observed properties of magne- 
tospheric accretion (Reiners, Basri & Christensen 2009). Re- 
cently, radio pulses were discovered from the L dwarf binary 
2MASSW J0746425+ 200032 with period P ^ 124 minutes, 
which point to a magnetic field of B ^ 1.7 kG (Berger et 
al. 2009) . The accretion rates in young BDs are smaller than 
in CTTS: M = 10~^^ - lO~^M0/yr, and are often strongly 
variable. For example, in 2MASSW J1207334-393254, M var- 
ied by a factor of 5 — 10 during a 6— week period (Scholz, 
Jayawardhana & Brandeker 2005) . We suggest that outflows 
may form in BDs during periods of enhanced accretion or in 
the propeller regime if the BD is rapidly rotating. 

As an example we consider a BD with mass Mbd = 
60Mj = O.O56M0, radius Rbd = O.IRq, and surface mag- 
netic field Bbd = 2 kG and obtain the reference parameters 
shown in the Table 1. The period of the star is another in- 
dependent parameter. Here we suggest rcor — 2 which cor- 



responds to P* = 0.13 days, which is a typical period for a 
BD. We also suggest that in Fig.[3]the star's radius is r = 0.5, 
that is, the disk is truncated at Vm = 1.2/0.5 = 2.4i?*. For 
these parameters we obtain an accretion rate of Mbd ~ 
1.8 X lO"^°M0/yr. For a smaller magnetic field, B = IkG, 
the same truncation radius will correspond to a smaller ac- 
cretion rate Mbd ~ 4.6 x lO^^^Mo/yr. The reference veloc- 
ity = 210 km/s is not different from the CTTSs case, and 
therefore the poloidal velocity of matter in the conical wind 
is Vp < (40 — 60) km/s. The higher- velocity component of the 
outflow, Vp ^ 200 km/s, can be easily explained if the BD is in 
the propeller regime. It is also possible that in the rotational 
equilibrium state the jet component is strong enough to drive 
jets. 

8.2 Application to Compact Stars 

8.2. 1 Symbiotic stars — white-dwarf hosting binaries 

Outflows are observed in some white-dwarf hosting systems. 
One class of them is the symbiotic stars (SSs). SSs are bi- 
nary stars in which a white dwarf orbits a red giant star 
and captures material from the wind of the red giant. Col- 
hmated outflows have been observed from more than 10 (out 
of ~ 200) symbiotic binaries. Most of them are transient and 
appear during or after an optical outburst that indicates an 
enhanced accretion rate (Sokoloski 2003). If SSs have a mag- 
netic field then enhanced accretion may drive conical-type 
outflows from the disk-magnetosphere boundary during peri- 
ods of enhanced accretion. The possibility of a magnetic field 
B ^6x lO^G in the SS Z And is discussed by Sokoloski & Bild- 
sten (1999) where flickering with a definite frequency was 
observed. In other SSs the magnetic field has not been esti- 
mated, but present observations do not rule it out (Sokoloski 
2003). The flickering in many SSs does not show a definite 
period, but the presence of a weak magnetic field is not ex- 
cluded (Sokoloski, Bildsten & Ho 2000). For a typical SSs ac- 
cretion rate ofM ^ lO^^M© /yr, a magnetic field as small as 
B ~ 3 X 10^ G will be dynamically important for the disk- 
star interaction. Thus it is possible that outflows are launched 
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from the vicinity of the SS as accretion-driven conical winds. 
ColHmation may be connected with a disk wind, disk mag- 
netic flux and/or the interstellar medium as discussed in §5.3. 



8.2.2 Circinus X-1 - the neutron-star hosting binary 

Circinus X-1 represents one of a few cases where a jet is seen 
from the vicinity of an accreting neutron star. The system is 
unusual because Type I X-ray bursts as well as twin-peak X- 
ray QPOs are observed. The neutron star is estimated to have 
a weak magnetic field (Boutloukos et al. 2006). The binary 
system has a high eccentricity (e ^ 0.4 — 0.9) and thus has 
periods of low and high accretion rates (e.g., Murdin 1980). 
Two-component outflows are observed. Radio observations 
show a non-stationary jet with a small opening angle on both 
arcminute and arcsecond scales. At the same time spectro- 
scopic observations in the optical (Jonker et al. 2007) and 
X-ray bands (laria et al. 2008) show that outflows have a con- 
ical structure with a half-opening angle of about 30°. Differ- 
ent explanations are possible for this conical structure, such 
as precession of a jet (laria et al. 2008). However, this ap- 
pears less likely because the axis of the jet has not changed 
in the last 10 years (Tudose et al. 2008). This neutron star 
may be a good candidate for conical winds, because (1) it has 
episodes of very low and very high accretion rates, and (2) 
a neutron star has only a weak magnetic field which can be 
strongly compressed by the disk, favouring the formation of 
conical winds. Table 1 shows possible parameters for neutron 
stars. Episodic collimated radio jets are also observed from the 
neutron-star hosting system Sco X-1 (Fomalont, Geldzahler, & 
Bradshaw2001). 



8.2.3 Application to black-hole hosting systems 

Jets and winds are observed from accreting black holes (BHs) 
including both stellar-mass BHs and BHs in galactic nuclei. 
The correlation between enhanced accretion rate and out- 
flows has been discussed extensively, and observational data 
are in favor of this correlation (e.g. Livio 1997). Recently, a 
conical-shaped ionized outflow was discovered in the black- 
hole hosting X-Ray Binary LMC X-1 (Cooke et al. 2008). It 
is not known what determines its shape, but the formation 
of conical winds is a possibility. Magnetic flux accumulation 
in the inner disk around the black hole was discussed by 
Lovelace et al. (1994) and Meier (2005) and observed in nu- 
merical simulations (Igumenshchev, Narayan, & Abramowicz 
2003; Igumenshchev 2008). Implosive accretion and outflows 
from black-hole hosting systems were analyzed by Lovelace et 
al. (1994) where angular momentum flows from the disk into 
a magnetic disk wind, leading to a global magnetic instability 
and strongly enhanced accretion. An accretion disk around 
a black hole may have an ordered magnetic field or loops 
threading the disk and corona. Fast accretion may lead to 
bunching of all field lines and possibly to conical winds. The 
inward advection of a large scale weak magnetic field thread- 
ing a turbulent disk is strongly enhanced because the sur- 
face layers of the disk are non-turbulent and highly conduct- 
ing (Bisnovatyi-Kogan & Lovelace 2007; Rothstein & Lovelace 
2008) . The mechanism of conical winds probably does not re- 
quire a special magnetic field configuration (such as a dipole) . 
Mohanty and Shu (2008) have shown that the X-wind model 
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Figure 23. Examples of HE I A10830 residual profiles and corre- 
sponding 1/2 veiling, r^, for 6 accreting CTTSs. The upper row shows 
high veiling objects with P Cygni profiles characterized by deep and 
broad blue absorption which is a sign of the high-velocity outflow 
The lower row shows low veiling objects with narrow blue absorp- 
tion which is a sign of a low-velocity outflow The latter also show 
red absorption from magnetospheric infall (from Edwards 2009) . 
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Figure 24. Modeling of the Hp line in RW Aurigae led to the conclu- 
sion that a cone-shaped wind with half-opening angle 30 — 40° and a 
narrow annulus gives the best match to the observations of this line 
(from Alencar et al. 2005). 



works when the star has a complex magnetic field configura- 
tion (see also Donati et al. 2006; Long, Romanova & Lovelace 
2007, 2008). 



9 CONCLUSIONS 

We have obtained long-lasting outflows of cold disk matter 
into a hot low-density corona from the disk-magnetosphere 
boundary in cases of slowly and rapidly rotating stars. The 
main results are the following: 

Slowly rotating stars (not in the propeller regime): 1. A 

new type of outflow — a conical wind — has been found and 
studied in our simulations. Matter flows out forming a conical 
wind which has the shape of a thin conical shell with a half- 
opening angle ^ 30°. The outflows appear in cases where 
the magnetic flux of the star is bunched up by the disk into 
an X-type configuration. We find that this occurs when the 
turbulent magnetic Prandtl number (the ratio of viscosity to 
diffusivity) Prm > 1, and when the viscosity is sufficiently 
high, av > 0.03. In earlier simulations of funnel accretion 
(e.g., R02; Romanova et al. 2003; Long et al. 2005) both vis- 
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Outburst: Approach to Equilibrium: Equilibrium: 

M^^^-lO-^M^yr' M _^ -10"^ M^^yr' M^^ -10"^ M^yr' 

Figure 25. Schematic model of an Exor VI 647 Ori. During the out- 
burst the accretion rate is enhanced, and the magnetospheric radius 
rm decreases and the magnetic field lines are bunched up (A). This 
results in a fast, hot outflow. As the accretion rate decreases, the 
disk moves outward and this results in a slower, cooler CO outflow 
(B) . Further decrease in the accretion rate leads to a quiescent state 
where the production of warm outflows stops (C). From Brittain et 
al. (2007) 

cosity and diffusivity were small and of the same order, and 
bunching of the magnetic field did not occur. 

2. The matter in the conical winds rotates with Keplerian ve- 
locity vk at the base of the wind and continues to rotate 
higher up. It gradually accelerates to poloidal velocities of 
Vp ^ 0.5vk. The conical wind is driven by the magnetic force 
which acts upwards and towards the axis. This is responsi- 
ble for the small opening angle of the cone, the narrow shell 
shape of the flow, and the gradual collimation of conical wind 
towards the axis inside the simulation region. 

3. Conical winds form around stars with different, including 
very low, rotation rates. The amount of matter flowing into 
the conical wind depends on a number of parameters, but 
in many cases it is Mw ~ (10 — 30)%Md. It increases with 
the rotation rate of the star and reaches almost 100% in the 
propeller regime. For rapidly rotating stars the outflows be- 
come strongly non-stationary. The period between outbursts 
increases with the spin of the star. 

4. There is another component of the outflow: a low-density 
high- velocity component of gas flowing along the stellar field 
lines. The volume occupied by this component increases with 
the rotation rate of the star. It occupies the entire region inte- 
rior to the conical wind in the propeller regime. 

5. A major part of the disk matter accretes to the star through 
the funnel flow and spins the star up. The conical winds carry 
away part of the disk angular momentum, but most of it is 
transported radially outward by the viscous stress. 

6. Conical winds can be further collimated at larger distances 
by the pressure of the surrounding medium (Lovelace et al. 
1991), disk winds (e.g. Fendt 2009), or by the magnetic flux 
threading the disk at large distances (e.g.. Matt et al. 2003). 

7. Conical winds may appear during strong enhancements of 
accretion, as in EXors or symbiotic variables. At the same time 
our simulations indicate that relatively steady outflows can 
exist for a long time (2 years in application to young stars) 
if the conditions for the magnetic field bunching are main- 
tained. 

8. Exploratory 3D simulations of conical winds from accret- 
ing stars with a significantly misaligned dipole field show that 
the conical winds are approximately symmetric about the ro- 
tational axis of the star (and the disk) . 



Propeller-driven outflows appear around rapidly rotating 
stars for conditions where rm > rcor and where the condi- 
tion for bunching, Prm > 1, is satisfied. Their properties are 
the following: 

1. Two distinct outflow components are found in the propeller 
regime: (1) a relatively low- velocity conical wind and (2) a 
high-velocity axial jet. 

2. A significant part of the disk matter and angular momen- 
tum flows into the conical winds. At the same time a signif- 
icant part of the rotational energy of the star flows into the 
magnetically-dominated axial jet. Formation of powerful jets 
is expected. This regime is particularly relevant to protostars, 
where the star rotates rapidly and has a high accretion rate. 

3. The star spins down rapidly due to the angular momentum 
flow into the axial jet along the field lines connecting the star 
and the corona. For typical parameters a protostar spins down 
in 3 X 10^ years. The axial jet is powered by the spin-down of 
the star rather than by disk accretion. 

4. The matter fluxes into both components (wind and jet) 
strongly oscillate due to events of inflation and reconnection. 
Most powerful outbursts occur every 1 — 2 months. The in- 
terval between outbursts is expected to be longer for smaller 
diffusivities in the disk. Outbursts are accompanied by higher 
outflow velocities and stronger self-collimation of both com- 
ponents. Such outbursts may explain the ejection of knots in 
some CTTSs every few months. Enhanced accretion due to ex- 
ternal factors will also lead to formation of a new blob/knot 
in the jet. 

The values of the transport coefficients ay and ad in real- 
istic accretion disks remain uncertain, but it is widely thought 
that they are due to Magneto-Rotational Instability (MRI)- 
driven turbulence (Balbus & Hawley 1998). MRI simulations 
suggest that the turbulence may give values of a in the range: 
= 10-2 _ 0.4 (e.g.. Stone et al. 2000). 

If in actual accretion disks the transport coefficients are 
large, a ^ 0.1, then strong outflows are expected during pe- 
riods in which ay > ad, or when the accretion rate is en- 
hanced. The condition for magnetic field bunching is that the 
magnetic Prandtl number of the turbulence Pr= ay /ad > 1. 
The effective Prandtl number may be significantly increased 
owing to the highly conducting surface layer of the disk 
(Bisnovatyi-Kogan & Lovelace 2007; Rothstein & Lovelace 
2008). The field lines may be bunched up not only due to 
high viscosity, but also by many other processes, e.g., due to 
thermal or magnetic instabilities, or due to the Rossby wave 
mechanism (Lovelace et al. 1999; Li et al. 2000). If the trans- 
port coefficients are very small, say ay ^ ad = 0.01, then 
quasi-stationary accretion through a funnel flow (with no out- 
flows) is expected (R02; Long et al. 2005). If the star rotates 
rapidly, it will be in the weak propeller regime in which a star 
spins down, but no outflows are produced, Romanova et al. 
2004a). If at some point the accretion rate is enhanced due 
to one or another mechanism, then conical winds will form in 
spite of the diffusivity being small. 
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APPENDIX A: VISCOSITY AND DIFFUSIVITY 

Here we discuss the treatment of viscosity and diffusivity in 
the axisymmetric code in greater detail. The stress tensor 
Tik —Tik + rik consists of an ideal part, 

-lik = pviVk -\-pOik -\- -—Oik z — , 

STT 4:7T 

and a viscous part nk which takes into account small-scale 
turbulent velocity and magnetic field fluctuations. 

We assume that the stress due to the turbulent fluctua- 
tions can be represented in the same way as the collisional 
viscosity by substitution of the turbulent viscosity coefficient 
into it. Moreover, we consider that the viscous stress is deter- 
mined mainly by the gradient of the angular velocity because 
the azimuthal velocity is the dominant component in the disk. 
The dominant components of the tensor nk in spherical coor- 
dinates are: 



-utpT sm 



dr 



dQ 
dO 
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Here — v^/r sin is angular velocity of the plasma and vt 
is the kinematic turbulent viscosity. 

Separating out the viscous stress in the component of 
eqn. ([2| gives 

d{pv^) , 1 d(r''Trcp) , 1 disin'OTe^) _ 



The azimuthal component of the induction equation gives 



dt 



+ 



1 d 



4 . .d^ 



dr 



+ 



r sm 

J ^ 

sin^ 6* dO 



p sm 



a6> 



dQ 



(Al) 



where T^^ and T^^ are components of the in viscid part of the 
stress tensor. 

The viscosity leads to dissipation of the kinetic energy 
and its conversion into thermal energy and to a correspond- 
ing increase of the entropy. In both types of runs (propeller 
and conical winds) we have neglected viscous heating. We 
have also neglected radiative cooling. Inclusion of heating 
and coohng is a separate and complex physics problem which 
is different for different types of stars. We suggest that the 
viscous heating is compensated by radiative cooling. Thus, 
the main "role" of viscous terms is the transport of angular 
momentum outward which allows matter to accrete inward 
to the disk-magnetosphere boundary. 

We also assume that the plasma has a finite magnetic dif- 
fusivity. That is, the matter may diffuse across the field lines. 
We assume that the finite diffusivity of the plasma is also due 
to the small-scale turbulent fluctuations of the velocity and 
the magnetic field. The induction equation averaged over the 
small-scale fluctuations has the form 



dt 



V X (v X B) + cV X = , 



(A2) 



Here, v and B are the averaged velocity and magnetic fields, 
and = — (v^ x B^) /c is electromotive force connected 
with the fluctuating fields. Because the turbulent electromo- 
tive force E^ is connected with the small-scale fluctuations, it 
is reasonable to suppose that it has a simple relation to the 
ordered magnetic field B. If we neglect the magnetic dynamo 
a-effect (Moffat 1978), then (v' x B') = -rytV x B, where 
r\t is the coefficient of turbulent magnetic diffusivity. Equation 
(|A2|) now takes the form 



aB 



V X (v X B) + V X (rytV x B) = . 



(A3) 



We should note that the term for E^ formally coincides with 
Ohm's law 



J= -^VxB = -^E^ . 
47r 47r?7t 

The coefficient of turbulent electric conductivity a — c? I^T^r\t. 
The rate of dissipation of magnetic energy per unit volume is 

£! = |i(VxBf . 
cr 47r 

To calculate the evolution of the poloidal magnetic field it 
is useful to calculate the ^-component of the vector-potential 
A. Owing to the assumed axisymmetry 

1 a(sin(9A^) ^ _ \d(rA^) 



B 



(A4) 



rsin^ dO ^ r dr 

Substituting B = V x A into the induction equation gives the 
equation for the component of the vector-potential 

dA^ fld^irA^) 1 d 1 d{smOA^ 



dt 



+ 



r'^dOsinO dO 



[vxB]^.(A5) 



dBd, 



r dr \ * dr 



1 d 



r]t d(sinOB^) 



dt r dr \ dr J dO \ sin 6 
1 (d[r{-v X B)g] a(v X B) 



r y dr dO 
The Joule heating rate per unit volume is 

4:7V 47vr^ 



(A6) 



(A7) 



djrBe) _ dBr 
dr dO 



+ 



djrB^ 
dr 



+ 



disinOB^) 



We included the Joule heating only in the propeller case runs, 
for completeness. However, we observed that although this 
term led to some heating, it was not the reason for the pro- 
duction the outflows in the propeller regime. Test simulations 
with no heating led to similar outflows, because the main 
driving forces are magneto-centrifugal forces. In the conical 
wind simulations we did not include Joule heating. Therefore 
we suggest that both Joule and viscous heating are exactly 
compensated by the radiative cooling. 



APPENDIX B: NUMERICAL METHOD 

For numerical integration of the MHD equations including 
the magnetic diffusivity and viscosity in the disk, we used a 
method of splitting of the different physical processes. Our 
simulation algorithm has a number of blocks: (1) an "ideal 
MHD" block in which we calculate the dynamics of the plasma 
and magnetic field with dissipative processes switched off; 
blocks (2) and (3) for the diffusion of the poloidal and 
azimuthal components of the magnetic field calculated for 
frozen values of the plasma velocity and thermodynamic pa- 
rameters (density and pressure); and block (4) for the calcu- 
lation of viscous dissipation in which we took into account 
only the rcj) and Ocj) components of the viscous stress tensor. 

(1) In the hydrodynamic block, the ideal MHD equa- 
tions are integrated numerically using an explicit conserva- 
tive Godunov-type numerical scheme. In our numerical code 
the dynamical variables are determined in the cells, while the 
vector-potential of the magnetic field, A^, is determined on 
the corners. For calculation of fluxes between the cells we use 
an approximate solution of the Riemann problem analogous 
to the one described by Brio & Wu (1988). 

For better spatial resolution, the restricted antidiffusion 
terms based on the MINMOD limiter are added to the fluxes 
(Kulikovskii, Pogorelov & Semenov 2001). The spacial split- 
ting has not been performed. Integration of the equations 
with time is performed with a two-step Runge-Kutta method. 
To guarantee the absence of magnetic charge, we calculate 
at each time-step the ^-component of the vector-potential 
Ad), which is then used to obtain the poloidal components of 
the magnetic field (Br, Be) in a divergence-free form (Toth 
2000). In other words, the condition V • B = is satisfied 
with machine accuracy. 

(2) In the block where the diffusion of the poloidal mag- 
netic field is calculated, we numerically integrate equation 
( [A5| ) for the ^-component of vector-potential. During this cal- 
culation we freeze the values of Ach on the inner and outer 
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boundaries of the simulation region. In the equatorial plane 
we have the symmetry conditions dAcf^/dO = 0. On the sym- 



metry axis we have = 0. Equation (A5) is approximated 
by an implicit difference scheme. The approximation is cho- 
sen so that the operator on the implicit time-layer is symmet- 
ric and positive. For solving the system of equations on the 
impHcit time-layer, we used ICCG (Incomplete Cholesky Con- 
jugate Gradient) method (Kershaw 1978). Because the size of 
the grid cells and the coefficient of magnetic diffusivity vary 
strongly in space, the elements of the matrix of the system 
also vary strongly. To remove this undesirable property, we 
changed the matrix so that it has diagonal elements equal to 
unity. 

(3) In the block where the diffusion of the azimuthal com- 
ponent of the magnetic field is calculated, we numerically in- 
tegrate equation ( [A6| ) . At the inner and outer boundaries, 

is frozen in this computational block. Along the rotation axis 
and on the equatorial plane = 0. Equation (A6 ) was ap- 
proximated by a numerical scheme with a symmetric positive 
operator on the implicit time layer. The corresponding system 
of linear equations is solved by the ICCG method. 

(4) In the block where the viscous stress is calculated, we 
integrate equation ( [Al| ) numerically for the angular velocity 
of matter Q = Vci)/r sin 0. At the inner boundary of the simu- 
lation region we take Q = Q*, the angular rotation velocity 
of the star. At the outer boundary, Q is taken to be fixed and 
equal to the corresponding Keplerian value. On the axis and 
in the equatorial plane we have the condition of zero stress for 
the Ocj) -component of the viscous stress tensor. Equation ( |A1| ) 
is approximated by a numerical scheme with a symmetric pos- 
itive operator on the implicit time-layer. The corresponding 
system of linear equations is solved by the ICCG method. 

The code has passed all the standard tests. In addition 
it has been used for the solution of a number of impor- 
tant astrophysical MHD flow problems (e.g., Ustyugova et al. 
1999; R02; R05; U06). In one problem, stationary super-fast- 
magnetosonic MHD outflows were obtained for the first time 
with the disk treated as a boundary condition (Ustyugova et 
al. 1999). In that work all of the flux function integrals of 
motion were calculated from the simulations and found to 
be constant as required by the theory (e.g., Lovelace et al. 
1986). Additionally, the cross-field force balance was checked 
numerically. Subsequently, similar results were obtained by 
Krasnopolsky et al. (1999) using the ZEUS code. This test of 
the simulations against the axisymmetric MHD theory is an 
important test in that it involves all three components of the 
flow velocity and the magnetic field. Simulation results of ac- 
cretion to a star with a dipole field (R02; Long et al. 2005) 
were recently confirmed by Bessolaz et al. (2008) . 



APPENDIX C: TESTS OF THE CODE 

We performed multiple tests of the code with different grids. 
We plan to describe the whole set of tests in a separate pa- 
per. However, in Appendix B we show two tests relevant to 
the ideal and diffusive blocks of the code. Below we show 
two examples of such tests. In the first test we checked the 
ideal MHD block of the code (with viscosity and diffusivity 
switched off). In the second test we checked the diffusion 
block of the code separately. 



CI Test of the ideal MHD block 

To check the MHD block of the code we performed the stan- 
dard "rotor problem" test. This test has been used by a num- 
ber of authors for testing MHD solvers including the energy 
equation (e.g., Balsara & Spicer 1999). We use this test to 
check our ideal MHD block of the code (viscosity and diffusiv- 
ity are switched off) with an isentropic equation dS/dt — ^ in- 
stead of the energy equation. In this situation the shock waves 
do not have physical sense. However, the goal of this test is 
not to test the physics of shock waves, but instead to demon- 
strate the ability of the numerical algorithm to solve the 2D 
adiabatic MHD equations. 

We solve the MHD equations numerically in the region 
—0.5 < X < 0.5, —0.5 < 2/ < 0.5 in a Cartesian geometry. At 
the beginning of simulations, t = 0, the pressure in the region 
is constant, p — I, and the magnetic field is homogeneous, 
Bx — 0, By — In the center there is a circle with radius 
ro = 0.1 (radius r — + y^^ where the density of matter 
is />o = 10 and the matter rotates as a solid body with angular 
velocity cjo = 20. At r > ri = 0.115, the density is pi = 1 and 
the matter is at rest. In the ring ro < r < ri, the density and 
velocity are linearly interpolated between those at r = ro and 
r — r\. 

The equations of ideal adiabatic MHD were solved with 
the Godunov-type scheme used in the "main" code. In the test 
we used a homogeneous grid with step Ax = A?/ = 1/A^, 
where = 100,200,400. The time-step is chosen from the 
condition At = 0.4Ax/^;macc, where Vmax is the maximum 
velocity of propagation of the perturbations. The results of 
the simulations are shown in Fig. |C1[ One can see that the 
density and the field line distribution is very similar in all 
three cases, while simulations at the highest grid resolutions 
give almost identical results and the convergence of the re- 
sults is evident. The bottom panels of Fig. [cT] show selected 
streamlines with numbers which confirm the similarity and 
convergence of the results. The test had been performed on a 
homogeneous grid, while our simulations were done in spher- 
ical coordinates with a high grid resolution near the star, and 
much coarser resolution at the outer boundary of the simu- 
lation region. In this paper we mainly investigate the launch- 
ing of jets and winds from the disk, and therefore we need 
to have adequate grid resolution in the region where mat- 
ter is launched from the disk into the winds, that is, at radii 
r ^ 2 — 3 for the conical winds, and r 3 — 5 for the propeller- 
driven winds. The grid resolution at these radii corresponds 
approximately to the homogeneous grid with the lowest grid 
resolution, 100 x 100. The above test shows that the grid res- 
olution in this region is sufficiently good for investigation of 
the physics of outflows in this region. It is clear that the grid 
resolution at larger distances, and in particular close to the 
outer boundary is not very high. However, this region does 
not influence the physics of the process. In this region mat- 
ter flows either inward due to the viscosity (in the disk), or 
outward as winds (in the corona). We expect larger simula- 
tion errors in these regions, but these do not change the main 
result: the launching of winds from the disk-magnetosphere 
boundary. The magnetic force accelerating matter into the jet 
component has higher accuracy close to the star and less ac- 
curacy at larger distances. However, Fig. lid shows that the 
main acceleration into the wind occurs close to the star where 
the grid resolution is good. 
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Grids 


max AA^ 


max AB^ 


51 X 31 


1.32 X 10-3 


7.11 X 10-4 


75 X 45 


6.20 X 10-4 


3.24 X 10-4 


101 X 61 


3.40 X 10-4 


1.75 X 10-4 



Table CI. The maximum absolute error in calculations of the and 
B(j) obtained in the diffusivity tests at different grid resolutions. 



We are planning future simulations with higher resolu- 
tions and in larger regions with the recently MPI-parallelized 
version of our code. 



C2 Test of the diffusivity block 

Here we test the diffusivity block of the code. For this we 
"switch off" the hydrodynamic fluxes and the correspond- 
ing right-hand side terms in equations A5 and A6 and inte- 
grate these equations numerically in the region ro < r < ri, 
< < 7r/2. In the case of the constant diffusivity coef- 
ficient (we take 77 = 1) equations A5 and A6 coincide, and 
so the problems for the azimuthal component of the vector- 
potential and for the azimuthal magnetic field differ only in 
the boundary conditions at the equator (at = 7r/2) and the 
initial conditions. That is, 



OA, 



0, 



0. 



0=7^/2 



6'=7r/2 

On the symmetry axis {_0 = 0) we have: 

A(f)\e=o — 0, B(f)\e=o — 0. 



We can find some particular solutions of equations (A5) and 
(A6) by the separation of variables method. For testing we 
choose the following solutions: 



C sin^ 



1 



5 sin 20 



1 + 



10ry(to - t) 



C sin cos ^ I _^ ^Tj(tQ — t) 



(CI) 



(C2) 



where C and to are constants. Equations for the 0- 
components of the vector-potential and magnetic field were 
integrated numerically using our 2.5D code (with the 
"switched-off" hydro fluxes and zero right-hand-side terms). 
We used the same grid geometry as in the main simulations; 
that is, spherical coordinates with an equidistant grid in the 
meridional direction and expanding grid in the radial direc- 
tion, where the expansion factor is determined by the fact that 
we keep all sides of each grid cell to be approximately equal. 
For the test we used the grid corresponding to the main sim- 
ulation runs in the conical wind, A^^ = 51 and Nq — 31, and 
also we used finer grids, 75 x 45 and 101 x 61. We took the in- 
ner and outer radii of the simulation region to be the same as 
in the conical wind case, ro = 1, ri = 15.9. For these bound- 
aries we set and B^p to be equal to the analytical values 
determined by equations [Cl] and |C2[ T he initial conditions 
corresponded to the exact solutions ( [Cl|a nd [C2| ) at t = 0. 

The constants C and to in |Cl| and |C2| were chosen so that 
the solutions at the inner boundary are of the order of unity. 
We used C — 0.01, to = 50. We integrated the equations up 
to ti 2to 100. 

Fig. [C2] (top panels) shows the result of the diffusivity test 
for the vector-potential, A^. The solid line in the left panel 
shows the numerical solution, A'T^'^, which is obtained by in- 
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a, ( , V ^ 




Figure C2. Test of the diffusivity block of the code on the grid 51 x 31. The top panels show analysis of the vector-potential, A^. The left hand 
panel shows the simulated value A^^^ (solid lines) and the exact solution A^^"^* (dashed lines). The middle panel shows the absolute error 
AA(f) = A^"^ — A^^^^^ , and the right hand panel shows the relative error, AA^/A^^"^*. The bottom panels show similar analysis but for the 
azimuthal component of the magnetic field, . 



tegrating equation A5, at time t = 100. The dashed line shows 
the exact solution, A^^"*""*, obtained by integrating the exact 
solution |C1| One can see that the solutions are very close to 
each other, and the dashed line is barely visible. The mid- 
dle panel shows the absolute error AA^ = A^^"^ - A^^"''*. 
One can see that the error is very small everywhere, taking 
into account the fact that the maximum value of the func- 
tion, A^ ^ 1. The error increases towards the star. However, 
this is because the value of the function A^ strongly increases 
towards the star. The right hand panel shows that the rel- 
ative error, (A^"""^ - A^^"''*)/A^^"''*, decreases towards the 
star. One can see that the relative error in calculation of the 
diffusivity is also small and is of the order of (1 — 2)%. Fig. 
C2 (bottom panels) shows similar analysis for the compo- 
nent, where equations (A6) and [C2] have been solved for the 
numerical and exact solutions. 

We performed similar simulations and analysis at finer 
grid resolutions, 75 x 45 and 101 x 61. Table 2 shows the 
maximum absolute error, max \A^'^ - A^^""""*!. One can see 
that at the finer grids, 75 x 45 and 101 x 61, this error is 2.1 
and 3.9 times smaller compared to the coarsest grid, 51 x 31. 
This confirms the convergence of the numerical solution to- 
wards the analytical solution. Similar comparisons for the 
component give factors 2.2 and 4.1, which also show conver- 
gence. 

A complete description of our methods and tests will be 
given in a separate paper. 



APPENDIX D: PARAMETER RANGES OF CONICAL WINDS 

To investigate the dependence on different parameters we 
took the main case and varied one parameter at a time. We 
performed several sets of runs: (1) with a fixed diffusivity and 
different viscosity coefficients; (2) with a fixed viscosity and 
different diffusivity coefficients; (3) with different rotation 
periods of the star P*; (4) with different coronal densities. 



Dl Dependence on viscosity at fixed diffusivity 

We fixed the diffusivity at = 0.1 and varied the viscosity 
coefficient in the range cxv — 0.01 — 1. Fig. |D1| shows matter 
fluxes onto the star and into the conical winds through the 
surface r = 6. We observed that for small viscosity, < 0.1, 
the magnetic field of the dipole diffuses through the inner re- 
gions of the incoming disk and an X-type configuration does 
not form. No conical winds appear in this case. We conclude 
that formation of conical winds requires ay > ad, that is, 
Prm > 1. Next we increased av and observed that an X-type 
configuration formed and conical winds were generated. We 
observed that the accretion rate to the star increases with ay, 
while the mass outflow rate into the conical winds increases 
but only slowly. For ay — 0.1, the matter fluxes onto the star 
and into the wind are small and approximately equal. For 
ay — 0.3 and ay — 0.4, the wind carries about 30% and 20% 
of mass respectively. The angular momentum carried to the 
star also strongly increases with ay . In all cases the star spins 
up, because the magnetospheric radius, rm ~ 1.2 is smaller 
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a =0.1 



a =0.1 




450 



Figure Dl. Matter and angular momentum fluxes onto the star and 
into the wind for fixed diffusivity ad = 0.1 but different coefficients 
of viscosity, av . 



a Star rotates more rapidly, the winds become more powerful, 
and the accretion rate to the star decreases, and therefore the 
situation becomes closer to the propeller regime. In rapidly 
rotating stars, the outbursts become episodic, and the interval 
between outbursts increases with spin. In addition, the vol- 
ume occupied by the fast coronal component increases with 
spin: at Vcor = 10 there is no fast component in the corona, 
at rcor = 3 it occupies some region above the conical winds 
(see Fig.[7]i), while at rcor = 2 it occupies a much larger re- 
gion (see Fig. [7|?) . At even higher spin this region occupies 
the whole simulation region as in the propeller regime. 
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Figure D2. Matter and angular momentum fluxes onto the star and 
into the wind for fixed viscosity av =0.3 but different coefficients of 
diffusivity, ad- 



than corotation radius rcor — 3. That is, the incoming matter 
brings positive angular momentum onto the star. The conical 
wind carries angular momentum away from the disk. Notice, 
however, that this is only a small part of the total angular 
momentum of the disk as shown above. 



D2 Dependence on diffusivity at fixed viscosity 

In the next set of runs we fixed the viscosity at = 0.3 
and varied the diffusivity: from 0.01 to 1. Fig. |D2| shows 
the integrated matter fluxes onto the star and into the wind 
at different diffusivities. We observed that no conical winds 
were formed for ad > dv. At relatively high diffusivity, 
ad — 0.1, 0.3, about 30% of the incoming matter flows into 
the conical wind. For = 0.1 (Prm = 3), the matter flux into 
the conical wind oscillates, while at ad — 0.3 (Prm = 1) no 
oscillations are observed. The conical wind also forms for very 
small diffusivity, ad — 0.01, but with slightly smaller matter 
flux into the wind. Angular momentum fluxes to the star and 
into the winds are approximately the same, excluding the case 
ad — 0.01 where the flux to the star is larger. 



D3 Variation of star's period 

We varied the period of the star (via the corotation radius) 
taking rcor — 3, 10 for slowly rotating stars and rcor — 
1.0, 1.5, 2 for more rapidly rotating stars. In the case of very 
slow rotation, rcor — 10, conical winds form and the out- 
flow rate into the wind is similar to that in the main case 
{rcor — 3), although the accretion rate onto the star is some- 
what larger. In stars with higher spin, rcor — 2, the amplitude 
of variability increases and the matter flux into the outflows 
increases up to 50% of the accretion rate to the star. For even 
higher spin, rcor = 1.5, the accretion rate to the star decreases 
by a factor of 5 compared with the main case, while the out- 
flow rate into the conical wind strongly increases. Thus when 



D4 Variation of coronal density 

Outflow of matter into a wind occurs if the corona is not very 
dense, and hence outflowing matter of the winds does not 
lose its energy while propagating through the corona. In the 
main simulation runs the initial density of the corona is 10^ 
times lower compared to the disk density { pc — 10 ~^ ver- 
sus pd — 10). To test the dependence on the coronal density 
we decreased its density by a factor of 3. These simulations 
showed that the matter fluxes onto the star and into the winds 
are not appreciably different from the main case. We conclude 
that the coronal density used in the main case is sufficiently 
small not to affect the outflows. 

The situation is different in the coronal region. For slowly 
rotating stars (the conical wind regime) there are no forces 
which tend to drive matter along the axis. Fast flow appears 
only in the part of corona where the stellar field lines are 
strongly inclined (see Fig. [s] Fig. [t]) . The rest of the corona 
has very slow motion towards or away from the star. The top 
half of the coronal region is matter-dominated (Fig. [Tsjz.d) 
and might be an obstacle to a fast outflow. We suggest that 
at lower coronal density the fast coronal component might 
occupy a larger region. 

In the propeller regime where the star rotates rapidly, the 
magnetic force is larger and drives the low-density coronal 
matter into the fast jet. The corona is magnetically-dominated 
(see Fig. [T6|i,d) during the whole simulation time. Initially, 
the density in the corona pc — 10~^ is 10 times lower than 
in case of slowly rotating stars. However, the initial density 
distribution is important only at the beginning of the long 
simulation runs. Later the density distribution is established 
by the outflow process. The main process is inflation of the 
dipole field lines. Inflating field lines carry matter along both 
disk and stellar field lines, and as a result, some matter pen- 
etrates into the corona. Strong disk oscillations and violent 
processes of inflation during which the conical wind compo- 
nent often changes its opening angle) lead to the penetration 
of a small amount of matter into the corona. The density in- 
creases away from the axis, where it is very small, towards the 
conical wind. In the case of a slowly rotating star, the disk os- 
cillations are weaker and there are no violent inflation events. 
For this reason the axial region above the star has a very low 
density. 



